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Preface 


The  purpose  of  this  study  was  to  assess  the  feasibility  of  detect¬ 
ing  turbulent  aircraft  wakes  uaing  Doppler  lidar.  Hardware  issues 
aside,  the  ultimate  detectability  of  an  aircraft  wake  is  determined  by 
the  level  of  backscatter  within  the  wake.  Therefore,  I  explored  the 
effects  of  three  different  scattering  sources  on  wake  detectability. 

The  three  scattering  sources  1  analyzed  were  ambient  atmospheric  aero¬ 
sols,  jet  engine  exhaust  soot  particles,  and  condensation  trail  (con¬ 
trail)  particles.  While  I  made  every  effort  to  be  as  complete  as 
possible  in  the  limited  amount  of  time  I  had,  by  no  means  is  this 
analysis  the  last  word.  In  particular,  I  found  that  enhancement  of  the 
backscattered  signal  by  the  presence  of  exhaust  soot  is  the  pivotal 
issue  on  whether  this  technique  will  work.  Unfortunately,  my  results  on 
soot  enhancement  are  only  a  rough,  first-cut  estimate.  A  good  deal  more 
research  is  required.  1  issue  this  statement  both  as  a  warning  to  the 
reader  and  as  a  suggestion  to  Interested  researchers. 

Of  course,  I  have  many  people  to  thank  for  their  help  on  this  proj¬ 
ect,  not  the  least  of  which  are  my  thesis  committee  members.  First  of 
all,  I  express  my  utmost  gratitude  to  Maj  (ret.)  James  Lange,  my  thesis 
advisor,  for  taking  me  on  and  for  taking  time  out  of  his  retirement  to 
assist  me.  In  addition,  I  extend  great  thanks  to  Dr  Won  Roh  for  his 
time.  Despite  the  fact  that  he  was  at  first  unfamiliar  with  coherent 
Ildar,  his  insightful  questions  and  suggestions  forced  me  to  think  out 
some  of  the  more  complex  issues  involved.  I  am  also  grateful  to  Dr  Ted 
Luke  for  taking  time  out  of  his  very  full  schedule  to  be  on  my  thesis 
committee . 

Many  more  people,  from  a  variety  of  organizations,  contributed 
their  time.  Of  these,  I  especially  wish  to  thank  Lt  Col  George  Koenig 
(AFGL/OPA)  and  Dr  Eric  Shettle  (Naval  Research  Lab)  for  their  helpful 
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comments  and  for  copies  of  BACKSCAT,  LOWTRAN,  and  the  MIE  scattering 
programs.  Their  help  was  instrumental  in  the  completion  of  this  thesis. 
I  also  extend  thanks  to  Mr  Ron  Rodney  (WRDC  Staff  Met)  for  running  the 
FASCODE  calculations.  Others  whose  contributions  are  appreciated  are  Mr 
Steve  Alejandro  (AFGL/OPA) ,  Dr  Joe  Wander  (HQ/AFESC),  Ms  Patrice  Acker¬ 
man  (Boeing),  and  Dr  Arnold  Barnes  (AFGL/LY) . 

Last,  but  far  from  least,  I  could  not  have  completed  this  thesis 
without  the  understanding  and  support  of  my  family.  Thankyou,  Barbl  I 
know  you  have  sacrificed  for  me.  And,  Kayla,  1  promise  I  will  have  more 
time  to  read  books  to  you  now  that  AFIT  is  almost  over. 

Michael  J .  Estes 
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Abie rice 


This  thesis  investigated  ths  feasibility  of  detecting  the  highly 
turbulent  air  within  high  altitude  aircraft  wake  vortices  using  infrared 
Doppler  lidar  mounted  in  a  search  aircraft.  This  technique  was  looked 
at  for  its  ability  to  detect  otherwise  stealthy  military  aircraft. 

Three  laser  wavelengths  were  analyzed:  1,064  /im,  2.091  /im,  and  9.113 
Mm.  Analysis  revealed  that  the  spectral  width  of  the  return  signal  from 
an  aircraft  wake  presented  a  good  signature  for  detection.  Based  or 
this  analysis,  a  minimum  .  lgnal-to-noisa  ratio  of  0  dB  was  established. 
Detection  performance  was  then  analyzed  using  signal- to-nolse  ratio  cal¬ 
culations  for  backscatter  by  ambient  atmospheric  aerosols,  jet  engine 
exhaust  soot  particles,  and  contrail  ice  particles.  Results  indicated 
that  atmospheric  aerosols  alone  were  not  sufficient  to  provide  enough 
backacattered  signal  for  detection  in  clean  regions  of  the  atmosphere. 
Backscatter  enhancement  by  exhaust  soot  particles  did,  however,  appear 
to  be  sufficient  for  detection  out  to  80  km  range.  Enhancement  by 
condensed  ice  particles  in  wake  contrails  provided  the  highest  signal 
levels  and  detection  well  beyond  100  km  in  range.  Interestingly,  the 
shorter  wavelength  Ildars  did  not  perform  nearly  as  well  as  the  9.115  /jm 
Ildar,  due  to  degradations  from  shot  noise,  wavefront  mismatch,  and 
refractive  turbulence. 
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DETECTION  OF  HIGH  ALTITUDE  AIRCRAFT  WAKE  VORTICES 
USING  INFRARED  DOPPLER  LIDAR:  AN  ASSESSMENT 

L _ Problem  Background 

Now,  more  chan  ever  before,  our  fighter  aircrafe  pilocs  muse  be 
keenly  aware  of  the  presence  of  approaching  enemy  aircraft.  To  this 
end,  long  range  detection  and  classification  of  possibly  hostile  air¬ 
craft  is  desired;  however,  recent  advances  in  stealth  technology  have 
made  aircraft  much  more  difficult  to  detect  by  conventional  means, 
namely  radar  and  passive  infrared  sensors.  At  the  same  time,  we  also 
desire  that  our  fighter  aircraft  become  more  stealthy.  That  means 
avoiding  the  use  of  radar,  which  immediately  flags  the  presence  and 
location  of  the  transmitting  aircraft.  In  order  to  maintain  control  of 
the  skies  in  future  conflicts,  we  need  to  explore  new  detection  schemes 
that  satisfy  both  of  these  requirements. 

One  possible  new  detection  scheme  that  shows  promise  is  the  detec¬ 
tion  of  aircraft  wake  vortices  using  Doppler  lidar.  As  a  by-product  of 
lift,  every  aircraft  in  flight  leaves  behind  a  turbulent  wake  of  air 
consisting  of  two  counter-rotating  vortices.  These  vortices  may  persist 
for  a  considerable  distance  behind  the  aircraft  before  breaking  up  and 
dissipating.  As  nn  example,  a  Boeing  707  aircraft  ir.  high  alititude 
cruise  leaves  behind  wake  vortices  that  persist  for  about  20  kilometers 
behind  the  aircraft  (Lutchen,  1984:13).  By  measuring  the  highly  turbu¬ 
lent  air  currents  in  the  wake,  it  may  be  possible  to  detect  the  presence 
of  an  approaching  aircraft  and  possibly  even  identify  it.  Doppler  lidar 
is  a  remote  sensing  instrument  that  appears  capable  of  detecting  such 
turbulent  airflow. 
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A  lidar  is  a  remote  sensing  instrument  used  to  probe  the  atmo¬ 
sphere.  Lidar  is  an  acronym  that  stands  for  Light  Detection  And  flang¬ 
ing  .  It  operates  much  the  same  way  as  conventional  radar,  except  that 
instead  of  transmitting  pulses  of  radio  waves  a  lidar  uses  a  laser  to 
transmit  pulses  of  light  waves.  The  transmitted  laser  pulse  scatters 
off  of  small  particles  (i.e.  -  dust  particles,  water  droplets,  mole¬ 
cules,  etc.)  suspended  in  the  atmosphere  and  a  small  fraction  of  this 
scattered  light  is  collected  by  the  lidar  receiver.  Doppler  lidar  is  a 
specialized  type  of  lidar  that  senses  the  Doppler  frequency  shift  in  the 
backscattered  light  caused  by  the  relative  motion  of  the  scattering  par¬ 
ticles  towards  or  away  from  the  lidar.  In  this  way,  a  Doppler  lidar  can 
measure  the  radial  air  flow  along  the  beam  path  and  thus  seems  to  be 
ideally  suited  for  detecting  aircraft  WAke  vortices. 

It  should  be  mentioned  here  that  lidar  may  also  be  used  to  detect 
trace  gas  species  in  tho  exhaust  trail  of  the  aircraft.  This  detection 
scheme  does  not  use  the  Doppler  shift  phenomenon  to  measure  air  flow, 
but  rather  it  uses  spectroscopic  lidar  techniques  such  as  DIAL  (Differ¬ 
ential  Absorption  Lidar),  Raman  scattering,  or  resonance  flourescence  to 
identify  specific  gases  peculiar  to  aircraft  exhaust.  Spectroscopic 
techniques  will  not  be  addressed  in  this  thesis;  however,  the  reader 
should  be  aware  of  their  existence  and  possible  usefulness  for  aircraft 
detection.  For  more  information  on  spectroscopic  lidar  techniques,  see 
Measures  (1984)  and  Hinkley  (1976) . 

The  concept  of  Doppler  lidar  detection  of  aircraft  vortices  has 
been  around  almost  as  long  as  the  laser  itself.  Possibly  the  first 
study  was  made  in  the  late  1960 's  and  employed  a  continuous  wave  C02 
laser-Doppler  velocimeter  to  detect  the  spinning  vortices  of  a  low  fly¬ 
ing  C-47  aircraft  (Huffaker  et  al,,  1970).  This  experimental  effort 
seemod  to  prove  the  viability  of  the  concept.  Since  that  time,  several 
more  studies  have  been  conducted,  some  of  which  were  concerned  with  low 
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altitude  vortex  detection  around  airports  for  the  purpose  of  aircraft 
safety.  While  related,  this  research  does  not  address  long  range  detec¬ 
tion  and  thus  is  not  of  great  interest  for  this  thesis.  Two  studies, 
however,  were  located  that  considered  long  range  detection  of  aircraft 
vortices  (Thomson  et  al.,  1989;  Lutchen,  1984).  One  problem  evident  in 
these  studies  is  that  maximum  detection  range  is  extremely  sensitive  to 
the  amount  and  type  of  atmospheric  aerosol  and  exhaust  particulate  pres¬ 
ent  both  in  the  aircraft  wake  Itself  as  well  as  along  the  laser  beam 
path.  The  purpose  of  this  thesis,  therefore,  is  to  explore  how  maximum 
detection  range  varies  with  atmospheric  conditions  and  aircraft  exhaust 
products . 

Three  infrared  laser  wavelengths  will  be  considered  in  this  analy¬ 
sis:  Nd:YAG  at  1.06  pm,  Ho:YAG  at  2.09  pm,  ar.d  isotopic  C02  at  9.11  /iin. 
Each  of  these  laser  systems  has  its  own  relative  merits.  The  Nd:YAG  and 
HoiYAG  lasers,  for  example,  are  both  solid  state  lasers  and  are  thus 
smaller,  simpler,  and  more  reliable  than  C02  gas  lasers;  however,  C02 
laser  technology  is  much  more  mature.  In  addition,  the  Nd:YAG  laser  is 
not  nearly  as  eyesafe  as  the  longer  wavelength  lasers.  These  relative 
advantages  and  disadvantages  will  not  be  addressed  in  this  thesis.  The 
only  criteria  to  be  compared  between  the  three  lasers  will  be  detection 
performance  as  a  function  of  wavelength. 

Maximum  detection  range  will  be  analyzed  for  backscatter  due  to 
naturally  occurring  atmospheric  aerosols  and  for  jet  engine  exhaust 
aerosols.  Each  of  these  aerosol  types  will  be  analyzed  separately. 
Atmospheric  aerosols  are  considered  here  to  be  non-cloud  particulates 
such  as  dust,  haze,  and  volcanic  debris  that  are  present  in  the  atmo¬ 
sphere  and  vary  with  altitude,  season,  geographic  location,  and  weather, 
The  type  and  amount  of  atmospheric  aerosol  affects  both  extinction  of 
light  along  the  beam  path  as  well  as  the  amount  of  light  backscactered 
to  the  lidar  receiver.  Exhaust  aerosols,  on  the  other  hand,  are  consid- 
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ered  to  be  by-products  of  combustion  and  thus  do  not  occur  naturally. 
Particles  in  the  exhaust  trails  of  aircraft  may  provide  enhancement  in 
the  backscattered  signal.  Exhaust  aerosols  to  be  analyzed  include  both 
soot  from  combustion  and  ice  particles  formed  in  condensation  trails 
(contrails).  The  magnitude  and  variation  in  all  these  effects,  which 
are  strongly  a  function  of  wavelength,  will  be  analyzed, 

Specific  signal  processing  algorithms  and  wake  physics  will  not  be 
covered  in  this  thesis.  While  overall  detection  performance  obviously 
depends  on  the  processing  algorithm  used,  the  ultimate  detectability  is 
limited  more  by  physical  constraints  such  as  backscatter  and  attenuation 
of  the  lidar  beam  as  well  as  the  sensitivity  of  the  Ildar  itself.  The 
analysis  of  detection  algorithms  will  be  left  for  future  research.  In 
addition,  wake  vortex  physics  is  a  complex  subject  and  will  only  be 
synopsized  here.  Only  simple  models  of  wake  vortex  behavior  will  be 
required  for  this  analysis. 

The  scenario  of  Interest  is  a  tactical  air-to-air  environment. 
Therefore,  the  most  probable  geometry  is  for  the  search  aircraft  (the 
one  with  the  lidar)  to  be  at  normal  cruise  altitude  of  from  20,000  to 
50,000  feet  with  the  lidar  looking  nearly  horizontally,  scanning  for 
aircraft  wakes  at  these  altitudes.  Ground-to-air  detection  and  air-to- 
air  detection  of  low  flying  aircraft  will  not  be  considered  here.  A 
study  by  Thomson  et  al.  (1989)  analyzes  air-to-air  detection  of  aircraft 
in  the  boundary  layer  and  lower  troposphere. 

This  thesis  is  organized  into  six  chapters.  The  first  is  the  cur¬ 
rent  background  chapter.  The  second  chapter  covers  theory,  including 
the  definitions  and  equations  to  be  used  in  succeeding  analyses.  The 
third  chapter  covers  preliminary  considerations  that  set  the  conditions 
and  assumptions  for  the  analyses.  The  fourth  and  fifth  chapters  present 
the  anaylsis  of  atmospheric  aerosols  and  exhaust  aerosols,  respectively. 


Each  aerosol  type  is  treated  Independently.  Finally,  the  sixth  chapter 
summarizes  the  results  end  discusses  recommendations  for  further 
research. 
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2_. _ Xheary 


This  chapter  presents  the  background  theory  to  be  used  in  subse¬ 
quent  analyses,  It  1b  assumed  that  the  reader  has  a  general  background 
In  optics  and  engineering,  but  is  perhaps  not  familiar  with  the  specific 
concepts  and  definitions  of  aircraft  wake  vortices,  atmospheric  optics, 
or  Doppler  Ildar.  The  first  subsection  in  this  chapter  deals  with  the 
general  structure  of  aircraft  wake  vortices.  The  second  subsection 
presents  a  quick  overview  of  applicable  atmospheric  optics  effects  such 
as  scattering  and  absorption.  Finally,  the  third  subsection  describes 
how  Doppler  Ildar  works  and  develops  equations  for  signAl-to-nolse  ratio 
and  velocity  estimation  error.  These  equations  will  bo  used  to  quantify 
detection  performance  and  thus  maximum  detection  range  for  the  three 
conditions  to  be  analyzed. 

IAj _ AlEgr.fl£.LJ,flkB  Vg.rtcx  Structure 

Every  fixed  wing  aircraft  in  flight'  generates  a  set  of  two  wake 
vortices  that  spin  off  of  the  wing  tips  and  form  in  parallel  behind  the 
aircraft.  Since  the  vortices  are  a  by-product  of  lift  by  the  wing,  they 
are  difficult  to  eliminate.  This  feature  makes  wake  vortices  good  tar¬ 
gets  for  detection  of  otherwise  stealthy  aircraft.  The  general  struc¬ 
ture  of  aircraft  wake  vortices  is  described  below  along  with  some  simple 
mathematical  models  of  airflow  in  the  vortices. 

As  an  aircraft  moves  through  the  air,  high  pressure  below  the  wing 
and  low  pressure  above  the  wing  generate  lift.  In  addition,  this  pres¬ 
sure  gradient  also  produces  air  circulation  about  the  wing,  which  even¬ 
tually  sheds  off  the  trailing  edge  of  the  wing  and  forms  into  two 
parallel  vortices.  This  process  illustrated  in  Figure  2-1.  Wake 
development  may  be  subdivided  into  three  stages:  1)  rollup  or  forma¬ 
tion,  2)  equilibrium,  and  3)  breakup  and  dissipation.  Figure  2-2  shows 
the  three  wake  stages,  which  are  discussed  further  below. 
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Figure  2-1,  Wake  Vortex  Formation  Behind  a  Lifting  Wing  (Lutchen, 

1984:6) 
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Figure  2-2.  Wake  Development  Stages  (Thomson  et  al.,  1989:3-2  as  taken 

from  Neilson  and  Schwind,  1971) 
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The  first  stage  of  wake  development  is  the  wake  roll  up  stage  dur¬ 
ing  which  the  trailing  vortices  start  to  form.  As  seen  in  Figure  2-1, 
the  circulating  "vortex  sheet”  rolls  up  behind  the  aircraft  forming  two 
counter-rotating  columns  of  air  parallel  to  the  line  of  flight.  Roll  up 
is  usually  complete  within  5  to  10  wingspans  behind  the  aircraft  and 
thus  is  not  of  much  interest  for  wake  detection  (Thomson  et  el., 
1989:3-1),  The  total  circulation  of  each  vortex  is  given  by 

r  (2-1) 

"pKtV. 

where 

T  -  circulation  parameter  (m*/s) 
hi  -  aircraft  mass  (kg) 

0  acceleration  of  gravity  (m/s*) 

P  -  air  density  (kg/m3) 
if  -  aircraft  wingspan  (m) 

K  -  wing  loading  constant 

-  it/  I  for  an  elliptically  loaded  wing 
V,  -  aircraft  velocity  (m/s) 

The  separation  distance  between  vortex  centers,  6,  is  proportional  to 
wingspan  by  the  wing  loading  constant,  or 

(>  -  A  s  ( 2  •  2 ) 

Once  the  wake  has  rolled  up,  it  stabilizes  into  a  relatively  time 
invariant  flow  structure.  This  stable  flow  field  is  the  second,  or 
equilibrium,  stage  of  wake  development.  The  length  of  the  equilibrium 
wake,  which  may  be  anywhere  from  100  to  1000  wingspans,  may  be  expressed 
as  (Thomson  et  al . ,  1989:3-9) 


where 

/.,  -  equilibrium  wake  length  (in) 
-  dimensionless  wake  length 


The  parameter  t,r  indicates  how  long  the  wake  remains  stable  before  ambi¬ 
ent  shear  and  turbulence  in  the  atmosphere  start  to  break  it  apart.  A 
model  of  the  dependence  of  on  atmospheric  conditions  may  be  found  in 
Crow  and  Bate  (1976).  For  this  analysis,  we  will  use  a  value  of  (^3 
(Thomson  et  al . ,  1989: 3-10). 

Several  models  exist  of  the  tangential  velocity  at  a  distance  r 
from  the  vortex  center.  A  model  by  Betz  (1932,  1933)  fits  experimental 
data  well,  but  has  a  discontinuity  at  the  vortex  center  where  velocity 
goes  to  Infinity,  Sprelter  and  Sacks  (1951)  proposed  a  model  to  elimi¬ 
nate  this  discontinuity.  The  Spreiter-Sacks  model  is  given  by  (Thomson 
et  al.  ,  1989:3-13) 


1,(0- 


/ 

\ 


Pur?' 

r 

kiw ' 


>■<>•, 

r  >  r. 


(2*4) 


where 

1,(0  -  tangential  velocity  component  (m/s) 
r  -  circulation  of  vortex  (m2/a) 
r  -  distance  from  vortex  center  (m) 

/ ,  -  vortex  core  radius  (in) 

A  plot  of  !',(/)  is  shown  in  Figure  2-3.  The  peak  tangential  velocity 


occurs  at  the  edge  of  the  center  region  of  each  vortex,  known  as  the 
vortex  core.  The  core  radius  is  given  by 


where  e^-0.197  (Spreiter  and  Sacks,  1951). 
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Radial  Distance  from  Vortex  Center  (m) 


Figure  2-3.  Wake  Vortex  Tangential  Velocity  Profile  for  a  Boeing  707 
Aircraft  Travelling  at  237  m/a  at  10.7  km  Altitude  (Lutchen,  1984:9) 
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The  flow  vector  at  any  point  within  the  equilibrium  wake  is  found 
by  the  vector  addition  of  the  flow  vectors  from  each  of  the  two  vor¬ 
tices.  A  cross  section  of  the  wake  with  contours  of  constant  velocity 
is  shown  in  Figure  2-4. 

The  elliptical  region  surrounding  both  vortices  (the  shaded  region 
in  Figure  2-4)  is  called  the  recirculation  cell.  The  dimensions  of  the 
this  region  are  1.64s  in  width  by  1.35s  in  height  (Lutchen,  1984:8). 

One  of  the  key  features  of  the  equilibrium  stage  is  the  fact  that  the 
size  of  this  cell  remains  fairly  constant,  The  recirculation  cell  is 
bounded  by  a  streamline  across  which  no  mass  is  transferred,  hence  the 
name.  The  fact  that  no  mass  is  transferred  across  this  sereninl  in<> 
implies  that  all  the  exhaust  emissions  become  entrained  in  the  recircu¬ 
lation  cell.  This  point  will  be  important  later  when  we  analyze  exhaust 
trails.  The  flow  at  the  edges  of  the  recirculation  cell  cause  the 
entire  cell  to  move  downward  with  a  vertical  sink  velocity  given  by 
(Lutchen,  1984:8) 

21  (2-6) 

In  addition  to  the  vertical  and  tangential  velocities  of  the  wake, 
axial  velocities  also  exist  due  to  thrust  and  drag.  The  vortex  cores 
have  a  forward  axial  velocity  due  to  drag  given  by  (Thomson  et  nl., 
1989:3-13) 

I  (2-7) 

1  - 

while  the  entire  recirculation  cell  has  a  backward  axial  velocity  due  to 
thrust,  which  is  given  by  (Thomson  et  al.,  1989:3-13) 

0.021’,  (2-8) 

1  1  "  (  1  .6-1 )  (  1  ,313) 
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SHADED  REGION  IS  RECIRCULATION  CELL 
*  -VALUE  OF  THE  STREAM  FUNCTION 


Figure  2-4,  Cross  Sectional  Flow  in  a  Stable  VHke  (Lutchen 
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where  both  I  ,  and  V,  are  in  units  of  m/s.  Again,  the  actual  vector- 

velocity  at  any  given  point  within  the  wake  is  the  vector  sum  of  all 
tangential,  vertical,  and  axial  velocity  vectors  at  that  point. 

Eventually,  interaction  between  the  vortices  and  turbulent  forces 
in  the  atmosphere  cause  the  wake  to  break  up  and  dissipate.  This  is  the 
third  and  final  wake  stage.  During  this  stage,  the  recirculation  cell 
begins  to  grow  and  dissipate  in  the  atmosphere,  According  to  the  liter¬ 
ature,  the  transverse  dimensions  of  the  wake  grow  at  a  time  rate  propor¬ 
tional  to  about  t1/3,  whore  t  is  time  (Thomson  et  al.,  1989:3-16).  Due 
to  its  unstable  nature,  we  will  not  be  concerned  with  the  third  wake 
stage  for  this  analysis, 

2.2. .  Atmospheric  Optics 

This  subsection  gives  a  quick  introduction  to  optical  scattering 
and  absorption  by  both  gases  and  aerosols  in  the  atmosphere.  Since  the 
subject  of  atmospheric  optics  is  very  complex,  only  the  basic  concepts 
are  presented  here.  In  addition,  we  will  only  consider  elastic  scatter¬ 
ing  processes  in  which  the  frequency  of  the  scattered  light  is  exactly 
the  same  as  the  frequency  of-  die  incident  light.  Wherever  possible, 
existing  models  of  these  atmospheric  processes  are  described  ,is  they 
will  be  used  extensively  in  the  analysis  sections. 

2.2.1.  ,  Preliminary  Definitions 

Consider  a  monochromatic  light  wave  incident  on  a  small  particle, 
portion  of  the  incident  light  will  be  scattered  off  in  all  directions 
(not  necessarily  equally),  a  portion  of  the  light  energy  will  be 
absorbed  by  the  particle,  and  a  portion  of  the  light  will  travel  on 
unaffected.  The  scattering  cross  section  of  the  particle  is  defined  as 
the  ratio  of  the  total  scattered  power  to  the  incident  irradiance  (Fenn 
et  al. ,  1985:18-15) ; 
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(2-9) 


where 

o,  -  scattering  cross  section  (m2) 

P,  -  total  scattered  power  (W) 

F0“  incident  irradiance  (W/m2) 

/,(0)  -  scattered  intensity  at  angle  O  (W/sr) 

Note  the  use  of  radiometric  symbols  £  for  irradiance  and  I  for  inten¬ 
sity.  Similarly,  for  absorption,  the  absorption  cross  section  of  the 
particle  may  be  defined  as  the  ratio  of  the  total  absorbed  power  to  the 
Incident  irradiance: 

l>„  (2-10) 


whore 

o„  -  absorption  cross  section  (m2) 

Pa  -  total  power  absorbed  (W) 

It  is  important  to  note  that  both  a,  and  o„  are  generally  strong  func¬ 


tions  of  wavelength,  as  are  all  the  parameters  that  will  be  defined  in 
this  subsection.  The  primary  consideration,  at  least  for  scattering,  is 
the  size  of  the  particle  relative  to  the  wavelength  of  light.  This 
wavelength  dependence  will  become  more  clear  in  the  following  sections 
on  molecular  scattering  and  aerosol  scattering, 

Now  consider  the  case  where  many  particles,  all  separated  with  ran¬ 
dom  spacing,  are  illuminated  by  the  light.  Assuming  low  particle  den¬ 
sity  and  single  scattering,  we  can  define  two  new  parameters  called  the 
scattering  and  absorption  coefficients  that  define  how  much  of  the 
incident  light  is  scattered  or  absorbed  per  unit  path  length  by  this 
group  of  particles  (Fenn  et  al.,  1985:18-15): 

(3,-A’o,  (2-11) 

-  A 'a„  (2-12) 
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where 


3,  -  scattering  coefficient  (m'1) 

P0  -  absorption  coefficient  (m"1) 

N  -  particle  density  (nf3) 

Furthermore,  we  can  define  the  extinction  coefficient  of  the  group  of 
particles  as  the  total  power  both  scattered  and  absorbed  per  unit 
length : 

P.-fl.  +  Pa  (2-13) 

Thus,  light  propagating  through  a  medium  of  particles  with  extinction 
coefficient  (3,  undergoes  an  incremental  decrease  in  irradiance  per  path 
length  dR  of 

dE  -  ~E[i,dR  (2-14) 

Integrating  this  equation  yields  the  Beer-Lambert  law  for  the  irradiance 
at  a  distance  from  the  source 

£(/?)- I(0)ex|j(-(i,/f)  (2-15) 

where 

f(/?)-  irradiance  at  distance  R  (W/m2) 

Note  that  the  equation  above  applies  only  to  a  uniform  medium  where  f3, 

is  constant  throughtout.  For  a  medium,  such  as  the  atmosphere,  where  [3, 
is  non-uniform,  the  Beer-Lambert  law  becomes 

(  r*  >  (2-16) 

/(ff)-f(O  pyp|  -  I  P, (<)((/'  1 
\  Jo  J 


The 


transmission  over  this  one-way  optical  path  of  distance  R 


is 

(2-17) 


For  a  single  scatterer,  the  intensity  of  the  scattered  light 
generally  depends  on  the  scattering  direction.  The  angular  scattering 
dependency  for  a  particle  Ls  summed  up  in  a  parameter  called  the  phase 
function,  which  defines  how  much  of  the  scattered  power  scatters  in  a 
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given  direction  and  in  a  given  solid  angle.  The  phase  function  is 
normalized  such  chat  integral  over  all  solid  angles  is  unity  (McCartney, 
1976:200) : 


-f 

4jVo 


/>(0)dcu  ■ 


(2-18) 


where 


/*(©)  -  normalized  phase  function  (sr*1) 

6  -  scattering  angle  (rad) 

For  a  volume  of  scatterers,  the  angular  scattering  coefficient 
indicates  how  much  light  will  be  scattered  in  direction  0  per  steradian. 
The  angular  scattering  coefficient.  P,(0),  is  the  normalized  phase  func¬ 
tion  per  particle  times  the  scattering  coefficient  for  the  ensemble  of 
particles.  In  symbols, 


(2-19) 


which  has  units  of  m^sr*1.  Note  that  the  previously  defined  quantity  (5, 

is  really  the  local  scattering  coefficient,  as  it  is  the  sum  of  the 
angular  scattering  coefficient  over  all  solid  angles 


(2-20) 


For  lidar  applications  we  are  interested  primarily  in  the  angular 
scattering  back  towards  the  light  source  at  0  -  ii .  The  angular 
scattering  coefficient  at  this  angle  is  also  known  as  the  volume 
backscatter  coefficient,  f3,( n )  and  is  given  as 


(2-21) 


Now  let  us  see  how  these  scattering  parameters  may  be  used  in 
actual  calculations,  For  a  single  scatterer,  the  scattered  intensity  at 
a  given  angle  may  be  expressed  in  terms  of  the  incident  irradlance  as 


2-12 


where 

F0  -  incident  irradiance  (W/m2) 

/,( 0)  -  acattered  intensity  <W/ar) 

For  a  volume  of  identical  scatterars  with  N  particles  per  unit  volume, 
the  scattered  intensity  is  given  by 

/.OO-P.O)*#!'  (2-23) 

•3,(9  )PtL 

where 

r  «■  volume  of  scatterers  illuminated  (m3) 

L  -  longitudinal  length  of  scattering  volume  (m) 

P o -  total  power  of  incident  light  (W) 

Take  note  of  the  fact  that  the  scattered  light  is  linear  in  terms  of 
power.  That  is,  the  scattered  power  or  intensity  is  a  linear  function 
of  the  number  of  particles.  This  observation  implies  that  the  scattered 
light  incoherently  adds,  even  though  we  have  made  no  mention  of  the 
coherence  properties  of  the  incident  radiation.  Indeed  this  is  the  case 
for  both  coherent  and  incoherent  illumination  if  we  assume  that  the 
number  of  scatterers  is  large  and  chat  they  are  randomly  spaced.  Since 
the  scattered  waves  from  all  particles  interfere  constructively  Just  at 
much  as  destructively  due  to  the  random  spacing,  the  interference  terms 
average  out  and  the  total  scattered  power  is  then  Just  the  sum  of  the 
power  scattered  from  each  individual  particle. 


WWWWhW-  Hsu  *  tl  »lWA\inWI*l  iM  *»  :  T  ■  v  f  WHH 


The  lower  80  kilometers  of  the  earth's  atmosphere  is  composed  of  a 
uniform  mixture  of  gases  consisting  of  79%  Nz,  21%  02,  and  much  smaller 
quantities  of  C02,  CHt ,  N20 ,  and  CO  (Fenn  et  al.,  1985:18*2).  Figures 
2-5, a-c  show  how  temperature,  pressuro,  and  density  vary  with  altitude. 
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Figure  2  5 


200 


rrpTr 

220 


IT  I  'H  I'TI  I  I 

240 


■I  ii  i  "|  n  nr 

260 


Temperature  (K) 


TrrrrrffTTn 

280 
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Figure  2-5, c.  Atmospheric  Density  versus  Altitude  (from  the  U.S.  Stan¬ 
dard  Atmosphere,  1976) 
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The  data  shown  in  these  figures  is  from  the  U.S,  Standard  Atmosphere 
(1976)  model.  The  scattering  of  light  by  atmospheric  gases  is  a  func¬ 
tion  of  the  index  of  refraction  of  air.  For  wavelengths  from  0.23  am 
through  the  infrared,  the  index  of  refraction  is  given  by  (Edlen,  1966) 

M7V1S.0)  (2-24) 


1  +  10  {  _°°*  1 -(v/b,)J+ 1 -(v/6,)’ J 


where 

n  -  index  of  refraction 
P  -  total  air  pressure  in  millibars  (mb) 
r o  -  1013.25  mb 

Pw  -  partial  pressure  of  water  vapor  (mb) 
7  -  temperature  (K) 

T0  -  273.15  K 

v  -  optical  wavenumber  (cm*1) 
ci  o  "  83  .  42 
ci,  -  185.08 
a t  -  4.11 


b  i  -  1 . 140xl05  cm'1 
b2  -  6.24x10*  cm'1 


c o  -  43.49 
c  ,  -  1 . 70x10*  cm*1 

The  dependence  of  optical  scattering  on  index  of  refraction  is  described 
by  Rayleigh  scattering  theory,  which  is  discussed  below. 

Rayleigh  scattering  theory,  developed  by  Lord  Rayleigh  in  1871  to 
explain  blue  sky  light  (McCartney,  1976:176),  describes  the  scattering 
of  electromagnetic  radiation  by  particles  whose  size  is  much  smaller 
than  the  wavelength  of  the  incident  radiation.  Such  is  the  case  of 
light  scattered  by  atmospheric  gas  molecules.  A  thorough  treatment  of 
Rayleigh  scattering  in  the  atmosphere  may  be  found  in  McCartney  (1976) 
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and  to  some  extent  in  Fenn  et  al .  ( 1985 : Section  18-1).  Here  we  will 
only  be  concerned  with  the  general  characteristics  of  Rayleigh  scatter¬ 
ing  and  how  to  compute  the  extinction  and  backscatter  coefficients. 

The  Rayleigh  scattering  coefficient  (total)  for  a  volume  of  gas 
molecules  is  given  as  (Fenn  et  al.,  1985:18*7) 

24n3{n*-  1  W6*3A\  (2-25) 

where 

\  -  wavlength  of  light  (m) 

N  •  particle  density  (m*s) 
n  -  refractive  index  of  the  gas 
A  -  depolarization  factor 

The  depolarization  factor  accounts  for  the  fact  that  air  molecules  are 
not  completely  isotropic.  The  depolarization  factor  for  dry  air  is 
given  as  A-0.0279  (Young,  1980), 

Note  that  0,  is  a  strong  function  of  wavelength.  One  of  the  main 


characteristics  of  Rayleigh  scattering  is  the  wavelength  dependence 
of  the  scattered  light  intensity.  Thus,  blue  light  scatters  much  more 
efficiently  than  red  light,  and  hence  the  blue  color  of  the  sky.  In 
addition,  we  see  that  the  amount  of  light  scattered  from  gas  molecules 
is  10,000  times  higher  at  a  wavelength  of  1  ^im  than  at  10  pm.  Other 
unique  characteristics  of  Rayleigh  scattering  (McCartney,  1976:176; 
include  the  fact  that  the  forward-  and  backscattered  intensities  are 
equal,  and  that  the  "light  scattered  at  90  degrees  is  almost  completely 
polarized" . 

The  angular  distribution  of  scattered  light,  or  phase  function,  for 
Rayleigh  scattering  is  given  by 


P(  0) 


3  2 

16)1(2*  A) 


IU 


A)*  (I  -  A)coo;'0] 


(2-26) 
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With  knowledge  of  the  scattering  coefficient  and  phase  function  we  may 
easily  compute  the  backscatter  coefficient  as 

(1,  (2-27) 

“  2K\'[ni*2j  \6-7Aj 

These  equations  thus  allow  us  to  calculate  the  angular  scattering 
of  a  volume  of  gas  of  index  of  refraction  n.  Absorption  of  light  by  gas 
molecules  is  a  much  more  complicated  subject,  as  seen  in  the  next 
section. 

2.2.1. _ Ab.a.o.r.p,Lifln. 

Like  molecular  scattering,  absorption  of  light  by  atmospheric  gases 
is  also  a  strong  function  of  wavelength,  although  in  a  much  more  compli- 
cated  way.  Each  gas  species  has  its  own  absorption  spectrum  consisting 
of  a  number  of  narrow  absorption  peaks  centered  about  the  electronic, 
vibrational,  and  rotational  energy  transitions  of  the  gas  molecules, 
Atmospheric  gas  absorption  is  then  the  combination  of  all  of  the  gas 
species  present  and  is  a  function  of  the  number  of  molecules  of  each 
gas,  the  wavelength,  and  the  temperature  of  the  air.  The  physics  to 
model  this  process  In  the  atmosphere  is  far  beyond  the  scope  of  this 
thesis;  however,  some  rather  accurate  computer  models  of  atmospheric 
absorption  have  been  developed,  most  notably  by  the  Air  Force's  Geophys¬ 
ics  Laboratory.  These  models  are  discussed  below  as  they  will  be  used 
in  later  analysis  sections. 

FASCODE  is  probably  the  most  extensive  model  of  atmospheric  absorp¬ 
tion  and  scattering.  FASCODE  takes  into  account  the  absorption  and 
scattering  characteristics  of  both  atmospheric  gases  and  aerosols  to 
calculate  transmission  and  radiance  along  an  optical  path.  FASCODE  com¬ 
putations  are  valid  in  the  spectral  region  from  the  microwave  to  the 
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near  ultraviolet  (Clough  et  al.,  1986:1)..  FASCODE  is  well  suited  for 
modeling  the  transmission  of  laser  light,  which  is  spectrally  very  nar¬ 
row.  FASCODE  will  be  used  to  compute  molecular  absorption  at  1.064, 

2  091,  and  9.115  pm  in  the  atmosphere. 

I+Ul. — Scattering  and. Absorption  by  Aerosols 

Aerosol  particles  suspended  in  the  atmosphere,  such  as  dust,  haze, 
smoke,  and  fog,  also  scatter  and  absorb  light.  These  particles  are  much 
larger  than  molecules,  typically  0.01  to  10  pm  in  diameter.  Scattering 
by  aerosol  particles  may  be  reasonably  modeled  by  Mie  scattering  theory. 
Mie  scattering  applies  to  particles  that  are  close  in  size  to  the  wave- 
longth  of  light.  Mie  scattering  theory  has  been  used  with  good  results 
to  modal  scattering  and  absorption  of  many  atmospheric  aerosols  (Fenn  et 
al. ,  1985:18-9  •  18-35) . 

Mie  scattering  theory  was  developed  around  the  turn  of  the  century 
by  Gustav  Mie  and  others  to  explain  the  scattering  and  absorption  of 
light  by  spherical  particles  of  arbitrary  index  of  refraction.  Mie 
theory  is  much  more  inclusive  than  the  Rayleigh  theory  in  that  it 
explains  scattering  by  any  size  particle  by  any  wavelength  of  electro¬ 
magnetic  radiation,  and  it  also  includes  absorption.  The  literature 
abounds  with  in-depth  treatments  of  Mie  scattering  theory,  which  is 
computationally  complex  (McCartney,  1976;  Blattner,  1972;  and  van  de 
Hulst,  1957).  The  following  is  a  quick  overview  of  Mie  scattering  that 
excludes  computational  details. 

In  Mie  theory,  the  absolute  values  of  wavelength  and  particle 
radius  are  unimportant.  Only  the  relative  ratio  of  particle  radius  to 
wavelength  matters.  Thus  a  convenient  parameter  called  the  size  parame¬ 
ter  is  defined  as 

,  ?.n  (2-28) 

a  *  k  r  ■  r 

K 
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where 


a  -  Mie  size  parameter 
Ac  -  2n/\  (m*1) 


r  -  particle  radius  (m) 

Consider  light  incident  upon  •  single  spherical  scatterer.  The 
scattered  radiation  field  consists  of  two  electric  field  components  with 
polarizations  parallel  and  perpendicular  to  the  plane  of  observation. 

The  angular  distribution  of  the  scattered  field  is  proportional  to  two 
dimensionless  field  distribution  functions,  St  and  S*,  which  are  in  turn 
functions  of  the  size  parameter,  a,  the  complex  index  of  refraction  of 
the  particle,  m,  and  the  scattering  angle,  G,  These  field  distribution 
functions  are  defined  as 


S,(a,  m ,  0)* 


I 

n  ■  I 


2n*  I 
n  ( n  4  ! ) 


(2-29) 


Sj(a  ,  m  ,  0)  * 
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2n*  1 
/i(/i  ♦  I ) 


(M. 


(2-30) 


where  the  functions  an.b„  are  the  complex  amplitude  coefficients 

computed  from  Ricatti-Bessel  functions  and  depend  on  the  size  parameter 
and  complex  index  of  refraction.  The  functions  Ji„ii'„are  computed  from 
Legendre  polynomials  and  are  only  dependent  on  the  scattering  angle 
(McCartney,  1976:229;  Blatcner,  1972:2).  The  field  distribution 
functions  are  related  to  two  intensity  distribution  functions  (also 
known  as  phase  functions),  which  are  proportional  to  the  scattered 
intensity  from  the  particle  at  a  givon  angle.  The  intensity 
distribution  functions,  denoted  by  i,  and  i8,  are  proportional  to  the 
scattered  intensities  of  the  polarization  components  perpendicular  and 
parallel  to  the  plane  of  observation  and  are  given  by 

l.-S.S'.-ISJ*  (2-31) 

<*-5js;-|S2|2  (2-32) 
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The  units  of  i,  and  t?.  are  sr'1.  At  scattering  angles  0-0  and  0«n,  i,  and 
<2  are  equal  (McCartney,  1976:264). 

The  total  intensity  backscattered  at  0*n  is  proportional  to  the 
incident  irradiance  as 


where 


4rt 


(2-33) 


f0"  incident  irradiance  (W/mJ) 

/,(»)  -  backscattered  intensity  (W/sr) 

It  may  be  seen  that  the  polarization  of  the  backscattered  light  is  the 
same  as  the  polarization  of  the  incident  light, 

The  total  scattering  cross  section,  ti„  and  extinction  cross  sec¬ 


tion,  o(,  of  tho  particle  are  given  by  (McCartney, 

h*0 


1976:247-248) 

(2-34) 


(2-35) 


The  p<  ''ceedlng  equations  assumed  a  single  scatterer  of  fixed  size. 
Now  consider  the  more  realistic  case  of  a  volume  of  particles  of  the 
same  index  of  refraction  but  different  sizes.  The  scattering  and 
extinction  coefficients  for  the  polydispersion  must  be  integrated  over 
the  size  distribution  of  particles  to  give  (McCartney,  1976:265) 

(2-36) 


o,(a)'K«)c/a 

i  r* 

P,  -  -  J  o,(a)/i(a)da 


(2-37) 


where 

n(a)  -  particle  size  distribution  (m3) 
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The  absorption  coefficient  is  again  defined  by 


Pn-P,-P, 

The  scattering  phase  functions  must  also  be  modified  to  incorporate  the 
integration  of  scattering  from  all  particle  sizes,  We  may  define  two 
new  phase  functions,  /* , ( Q )  and  P2C©).  for  the  polydispersion  of  scatters 
in  terms  of  the  phase  functions  for  each  individual  particle  as 
(McCartney,  1976:265) 

4 n  f°>  (2*38) 

n(a)i/(fl.a)da 

for  y-1,2 

As  before,  0)  and  /'f( 0)  are  proportional  to  the  scattered  intensity  at 


an  angle  O  per  unit  solid  angle  for  the  polarizations  perpendicular  and 
parallel  to  the  plane  of  observation,  Since  the  backscattered 
polarization  is  the  same  as  the  incident  polarization,  /*,(ii  )  -  /*2( n ),  and 
the  backscatter  coefficient  for  the  polydisparsion  of  particles  is  given 

by 


P.(«) 


Hi  /^c*1) 

4it[  2 


(2-39) 


As  mentioned  previously,  Mie  scattering  theory  applies  to  spherical 
particles.  In  general,  though,  atmospheric  aerosols  are  not  spherical. 
Liquid  aerosols,  such  as  water  droplets,  are  approximately  spherical, 
but  dry  aerosols  like  dust  and  soot  are  definitely  not.  Some 
researchers  feel  that  this  discrepancy  is  not  really  a  problem  since 
most  particle  sizing  instruments  tend  to  assign  an  "equivalent"  particle 
radius  for  a  spherical  particle  with  similar  scattering  properties  ns 
the  actual  particle,  The  resulting  models  of  atmospheric  aerosols 
derived  from  these  measurements  therefore  should  remain  approximately 
valid  (Fenn  et  al.,  198 5:18-15) .  In  the  specific,  case  of  cirrus  clouds, 
however,  the  spherical  approximation  gives  generally  poor  results 
(Shettle,  1990).  In  addition,  depolarization  from  cirrus  clouds  may  be 
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substantial.  That  is,  the  polarization  of  the  backseat tered  light  may 
be  significantly  different  from  the  polarization  of  the  incident  light. 

In  this  erse,  the  single  parameter  backscatter  coefficient  previously 
defined  does  not  apply  and  scattering  matrices  must  be  used  to  account 
for  polarization. 

2.J- _ Dflanlsi-.il.flar 

By  measuring  the  Doppler  frequency  shift  in  backscattered  light,  a 
Doppler  lidar  may  be  used  to  measure  the  relative  motion  of  air  cur¬ 
rents,  The  technique  most  commonly  used  for  this  measurement  is  hetero¬ 
dyne  detection,  which  is  described  below. 

A  complete  and  thorough  treatment  of  Doppler  lidar  does  not  really 
exist  in  the  open  literature  since  the  topic  pulls  in  expertise  from  a 
number  of  diverse  fields  such  as  radar,  laser  physics,  atmospheric  phys¬ 
ics,  etc,  Menzies  and  Hardesty  (1989),  though,  do  give  an  excellent 
overview  of  Doppler  lidar  for  atmospheric  wind  measurements  and  include 
a  number  of  classic  references  on  the  subject. 

2AL _ fiflD&r.al  PescrlRBlHD 

Most  Doppler  lidars  use  heterodyne  detection  to  obtain  the  very 
fine  spectral  resolution  required  to  measure  small  Doppler  shifts.  Such 
lidars  are  generally  monostatic  or  coaxial,  meaning  that  the  transmitted 
and  received  beams  use  the  same  optics.  Figure  2-6  shows  the  general 
layout  of  a  monostatic  heterodyne  lidar  system.  A  highly  coherent 
pulsed  laser  generates  the  output  beam,  which  is  expanded  by  a  telescope 
and  transmitted  out  into  the  atmosphere.  Backscattered  light  is  then 
collected  by  the  telescope,  photomixed  with  a  stable  reference  beam  from 
a  local  oscillator  laser,  and  then  focused  down  onto  a  square-law  photo¬ 
detector.  The  local  oscillator  frequency  is  slightly  offset  from  the 
transmitter  frequency  so  that  the  irradiance  on  the  detector  is 
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Figure  2-6.  General  Layout  of  a  Monostatic  Heterodyne  Lidar 
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modulated  by  the  difference  or  beat  frequency  between  the  reference  and 
received  beams.  The  detector  output  then  consists  of  a  beat  signal 
whose  mean  frequency  is  proportional  to  the  average  wind  velocity  within 
the  scattering  volume  and  whose  frequency  spread  is  proportional  to  the 
wind  variability  or  turbulence.  The  electronic  bandwidth  of  the  detec¬ 
tor  and  downstream  amplifiers  filtor  the  output  signal  to  a  finite  band¬ 
width  B. 

The  Doppler  frequency  shift  A/0  in  the  backscattered  signal  due  to 
a  wind  velocity  c  along  the  beam  path  is  given  by 


where  k  is  the  wavelength  of  the  transmitted  light,  The  resulting  bent 
frequency,  ,  is  then 

(2*41) 

where 

/ r# /  “  local  oscillator  optical  frequency  (Hz) 

/,  -  transmitter  optical  frequency  (Hz) 

Microscale  turbulenc :  within  the  scattering  volume  generates  a  distribu¬ 
tion  of  particle  velocities.  The  resulting  backscattered  signal  has  a 
Doppler  frequency  "spread"  due  to  this  distribution.  The  spectral  width 
of  the  received  signal,  then,  is  due  to  the  combined  spectral  widths  of 
the  transmitted  pulse  beam,  the  local  oscillator  beam,  and  the  backscat¬ 
tered  light.  The  spectral  widths  of  the  transmitter  and  local  oscilla¬ 
tor  beams  are  inversely  related  to  the  coherence  times  of  these  beams. 
Therefore,  in  order  to  measure  wind  velocity  with  any  accuracy,  the 
coherence  times  for  these  lasers  should  be  long.  This  requirement  loads 
to  the  commonly  used  term  coherent  Ildar  in  reference  to  heterodyne 
detection  lidar. 
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Since  the  transmitted  beam  is  pulsed,  range  information  may  be 
extracted  from  the  received  signal.  If  the  output  pulse  is  transmitted 

at  time  *-0,  then  the  received  signal  at  time  t  is  due  to  a  scattering 
volume  at  a  range  of  *c(,  where  c  is  the  speed  of  light.  The  scattering 
volume  has  transverse  dimensions  equal  to  the  beam  spot  size  at  that 
range  and  longitudinal  dimension  |cr,  where  t  is  the  transmitted  pulse 
length.  Thus  the  lidar  cannot  resolve  scattering  volumes  shorter  than 
one  half  the  pulse  length.  The  minimum  resolved  scattering  volume  is 
often  referred  to  as  a  range  gate  interval. 


In  this  subsection,  we  will  look  at  the  key  properties  of  the  back- 
scattered  light  that  is  incident  upon  the  lidar  receiver.  We  may 
describe  this  optical  signal  statistically  in  terms  of  its  power, 
frequency,  and  coherence. 

The  electric  field  of  the  backscattered  signal  is  n  sum  of  the 
electric  fields  from  a  large  number  of  independent  ecatterers.  Since 
the  scetterers  are  randomly  distributed  throughout  the  scattering  vol¬ 
ume,  the  resulting  electric  .field  is  a  random  phasor  sum  whose  amplitude 
is  Rayleigh  distributed  and  whose  phase  is  uniformly  distributed 
(Goodman,  1985:44-50),  In  terms  of  instantaneous  power,  the  return  sig¬ 
nal  is  then  an  exponentially  dlotributed  random  variable  with  equal  mean 
and  standard  deviation.  The  relative  uncertainty  in  power  is  therefore 
100%.  For  A'  independent  measurements ,  the  relative  uncertainty  falls 
to  100 \/{N  (Hardesty  et  al . ,  1981:3768-3769). 

The  average  optical  power  collected  by  the  lidar  receiver  from  a 
scattering  volume  at  range  R  may  be  found  using  radiometry.  Since  we 
are  primarily  interested  in  long  range  detection,  we  may  assume  that  the 
scattering  volume  is  in  the  far  field  of  the  lidar  receiver.  That  is, 
we  assume  that  the  scattering  volume  is  essentially  a  point  target.  If 
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the  transmitted  pulse  begins  at  t“ 0  and  ends  at  then  the  total 

power  collected  by  the  receiver  is  the  convolution  of  the  pulse  temporal 

power  profile  with  the  range  dependent  bAckscatter  coefficient  and 

tranamiasion  (Menzias  and  Hardesty,  1989:452).  In  symbols, 

r>f'  (  f  a  '  (2-42) 

M O-n.j  P,  (-2-  (!,(#)  -j  It \R)dR 

where 

n„  -  collection  efficiency  of  optics 

Pr(0  -  total  optical  power  incident  on  receiver  (W) 

P,(0  -  transmitted  pulse  temporal  power  profile  (W) 

-  volume  backscatter  coefficient  (m^sr"1) 

A  -  area  of  receiver  entrance  aperture  (m2) 
t(/f)  -  one-way  transmission  to  range  R 
R  -  range  to  scattering  volume  (m) 

Note  the  change  in  notation  on  the  backscatter  coefficient  from  previous 
definitions,  The  notation  used  here  is  more  common  in  the  lidar  commu¬ 
nity,  while  the  previous  notation  is  mors  common  in  the  optical  scatter¬ 
ing  community.  The  backscatter  coefficient  and  one-way  transmission 
have  been  written  us  functions  of  range  to  underscore  their  strong  range 
dependence . 

If  we  assume  that  P„(/?).  t (P),  and  R  2  do  not  change  appreciably 


over  the  range  gate  interval  between  •  c ( /  - t „ )  and  ^ r / ,  then  these  terms 
may  be  pulled  out  of  the  integral  above  leaving  only  the  pulse  temporal 
profile  term.  The  integral  of  the  pulse  power  over  the  pulse  length 
just  gives  the  total  energy  f,  in  the  pulse  times  one  half  the  speed  of 
light.  Note  that  f  here  denotes  energy,  not  irradiance  as  in  previous 
definitions.  The  simplified  expression  for  received  average  power  is 
then 


R,(0 


(2-43) 
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The  backscatter  coefficient  (3„(£)at  range  R  is  the  sum  oi  the 

backscatter  coefficients  for  molecules  and  aerosols,  or 

<2‘44> 

where 

“  molecular  backscatter  coefficient  (m*lar*1) 

PJ(/?)  -  aerosol  backscatter  coefficient  (m^sr'1) 

Likewise,  the  extinction  coefficient  (!,(#)  that  contributes  to  the  one* 

way  path  transmission  x(R)  at  range  R  is  also  the  sum  of  molecular  and 
aerosol  extinction  coefficients: 

P.<*)-Pr<*)*PJ(#>  (2‘U5) 

where 

p;‘(tf)  -  molecular  extinction  coefficient  (m*1) 

(?!'(/?)  -  aerosol  extinction  coefficient  (m* 1 ) 

We  noted  above  that  the  instantaneous  power  in  the  backscnttered 
signal  is  an  exponentially  distributed  random  variable.  As  the  par¬ 
ticles  move  in  random  motion  within  the  scattering  volume,  the  instanta¬ 
neous  power  changes.  The  rate  at  which  signal  power  fluctuates  is  a 
function  of  particle  velocities.  After  a  short  period  of  time,  the 
particles  have  rearranged  themselves  enough  that  the  new  instantaneous 
signal  power  is  independent  of  the  previous  power.  In  addition,  the 
phase  of  the  return  signal  is  independent  of  the  earlier  phase;  hence 
the  two  signals  are  uncorrelated  or  temporally  incoherent.  The  time 
period  over  which  the  return  signal  remains  coherent  (at  least  par¬ 
tially)  is  the  coherence  time  and  depends  on  the  distribution  of  par¬ 
ticle  velocities. 

Using  the  Wioner-Khinchin  Theorem,  the  power  spectrum  of  the  return 
signal  is  related  to  the  autocorrelation  of  the  signal  by  the  Fourier 
transform  (Goodman,  1985:74).  Thus,  high  particle  velocities  cause 
short  coherence  times  or  narrow  autocorrelations,  which  in  turn  mean 
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broad  power  spectrums.  Low  particle  velocities  correspond  to  narrow 
power  spec crums .  Assuming  perfectly  monochromatic ,  and  hence  infinitely 
continuous,  illumination,  the  shape  of  the  power  spectrum  is  the  same  as 
the  distribution  of  particle  velocities.  Therefore,  we  may  interpret 
the  spectrum  as  the  distribution  of  Doppler  shifte  from  the  volume  of 
particles. 

At  this  point,  let  us  consider  backscacter  from  molecules  and  back- 
scatter  from  aerosols  separately.  Host  atmospheric  aeroeols  are  small 
snough  that  thsy  movs  frssly  with  the  wind,  and  thus  the  mean  Doppler 
shifts  from  aerosols  end  molecules  within  e  scattering  volume  are  equiv¬ 
alent.  The  distributions  of  velocities,  however,  are  very  different. 

Molecular  velocity  distribution  is  a  function  of  absolute  tempera¬ 
ture  as  given  by  the  Maxwell  velocity  distribution  (Tlpler,  1987:69). 

The  velocity  distribution  in  one  direction,  say  the  x-direction,  is  a 
Gaussian  diatribution  with  standard  deviation 

[FT-  (2-46) 

u  V  Al.< 

where 

ou  -  velocity  etendard  deviation  (m/e) 

R'  -  universal  gaa  constant  (8.314  J/mole*K.) 

T  -  absolute  temperature  (K) 

M  ,  -  molar  mass  of  air  (kg/mole) 

McCartney  (1976:52)  gives  the  molar  mass  of  air  as  Afr28.964  g/mole. 

Figure  2-7  shows  a  plot  of  ou  versus  altitude.  As  the  plot  indicates, 
molecular  velocitiea  in  the  atmosphere  are  fairly  high.  Even  the  very 
cold  tropopaute  region  sees  molecular  velocities  of  *250  m/s.  The 
resulting  power  spectrum  from  moleculsr  backscatter  is  very  broad. 

Because  of  their  larger  size,  aerosols  have  a  typical  velocity 
spread  of  only  about  1-2  m/s  (Menzies  and  Hardesty,  1989:454),  and  the 
insulting  power  spectrum  is  relatively  narrow,  The  range  of  mean  wind  s 
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Altitude 


Relative  Power 


Figure  2-8.  Typical  Power  Spectrum  of  Molecular  and  Aerosol  Returns 
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peeds  in  the  atmosphere  is  usually  well  within  *'0  m/s .  Even  the  peak 
velocity  in  the  707  wake  vortex  is  only  18  m/s .  Over  the  range  of  mean 

air  velocities,  the  power  spectrum  of  molecular  backscatter  is  approxi¬ 
mately  flat.  For  this  reason,  molecular  backscatter  should  really  be 
considered  as  noise  power,  not  signal  power.  Signal  power,  therefore 
only  comes  from  aerosol  backscatter.  Figure  2*8  illustrates  the  rela¬ 
tionship  between  aerosol  and  molecular  Doppler  apectrums. 

The  spatial  coherence  of  the  backscattered  signal  is  also  degraded 
by  the  atmosphere.  Random  fluctuations  in  the  refractive  index  of  air 
between  the  lidar  and  the  scattering  volume  are  known  as  refractive 
turbulence.  The  refractive  index  structure  constant,  Cl,  indicates  the 


strength  of  refractive  turbulence  and  is  given  by 


(2-47) 


where 


it ,  -  refractive  index  at  point  1 
itj  -  refractive  index  at  point  2 
L  -  distance  between  points  1  and  2  (m) 


C?  has  units  of  m'a/3.  The  bur  over  the  numerator  in  the  previous  equa¬ 


tion  denotes  ensemble  average. 

The  effect  of  propagating  light  through  a  turbulent  medium  is  a 
reduction  in  the  spatial  coherence  of  the  wavefront.  The  transverse 
coherence  length,  /,,,  is  a  measure  of  spatial  coherence  as  it  indicates 
the  distance  along  the  wavefront  over  which  the  wavefront  is  coherent. 
The  transverse  coherence  length  of  the  backscattered  signal  is  given  by 
(Huf faker  and  Targ,  1988:66) 

.  T  r •  I'3'9  (2-48) 

r„  ■  0.058 1  C\0  )rs,iclr 

l  •'o 
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Transverse  Coherence  Length  (m) 


where  R  is  the  range  to  the  scattering  volume.  All  lengths  in  this 
equation  are  in  units  of  meters. 

At  nominal  aircraft  cruise  altitudes  between  5-15  kilometers  (16-50 
kft),  C\  is  fairly  constant  at  around  10'17  m'*7*  (Fenn  et  al., 

1985:18-72;  Goodman,  1985:429).  This  value  will  be  used  in  subsequent 
analyses.  Figure  2*9  shows  the  transverse  coherence  length  plotted  ver- 
sus  range  over  a  horizontal  path  for  the  three  Infrared  wavelengths  on 
interest , 

2.3.3. _ Signal- to-Nolaa Ratio 

Noise  limits  the  ultimate  sensitivity  of  any  lidar  system.  There¬ 
fore,  the  ratio  of  signal  power  to  noise  power  in  the  output  signal 
gives  a  good  one-parameter  indicator  of  the  performance  capability  of  a 
given  lidar  under  a  given  set  of  conditions.  Here  we  will  develop  an 
equation  to  express  the  signal- to-noise  ratio  of  a  heterodyne  detection 
Doppler  lidar. 

After  collection  by  the  telescope,  the  received  light  is  photomixed 
with  the  local  oscillator  reference  beam  and  focused  down  onto  the 
detector.  The  resulting  heterodyne  optical  signal  incident  on  the 
detector  has  the  form  (Post,  1985) 

+  (2'49) 

where 

/’,(<)»  optical  power  incident  on  detector  (W) 

-  local  oscillator  power  (W) 

P r  -  return  signal  power  (W) 

i|„  -  heterodyne  mixing  efficiency 

<xi, „  -  local  oscillator  optical  frequency  (rad/s) 

uu,  -  received  optical  frequency  (rad/s) 
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Following  the  previous  discussion  about  molecular  and  aerosol  signal 
characteristics,  only  the  aerosol  backscatcer  contributes  to  meaningful 
signal.  If  we  denote  the  optical  power  from  aerosols  as  ft,,  then  the 
signal  currant  t,(0  out  of  the  detector,  with  D.C.  components  filtered 
out,  has  the  form 


where 


r|„  -  detector  quantum  efficiency 
e  -  electron  charge  (C) 
h  -  Planck's  constant  <6 . 626 • 10-54  J  s) 
v  -  optical  frequency  (Hz) 

We  may  now  find  an  expression  for  the  noise  power  in  the  output 
current,  Shot  noise  in  the  output  signal  is  due  primarily  to  local 
oscillator  power  and  is  given  by 

<<!>«.,  "2o/li  (2-51) 


where 


<  «,*  -  average  shot  noise  power  (A*) 

i  -  average  detector  current  (A) 
ft  -  electronic  bandwidth  (Hr.) 

Shot  noise  usually  dominates  thermal  noise,  and  thus  thermal  noise  is 
neglected.  The  noise  contribution  from  molecular  backscatter  may  be 
considered  background  noise.  Background  noise  power  is  then 


<<!>*«*■ 


(2-52) 


where 

<  > «>ac*  “  average  background  noise  power  (A1) 

ft c  «  optical  power  from  molecular  backscatter  in  bandwidth  ft  (W) 
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The  total  noise  power  is  the  sum  of  all  noise  power  terms.  At  the 
wavelengths  of  interest,  however,  shot  noise  dominates  background  noise 
from  molecular  scattering.  Simulations  show  that  background  noise  is 
less  that  2%  of  shot  noise  for  a  lidar  located  at  9  kilometers  altitude 
looking  horizontally.  Thus  we  may  neglect  background  noise. 

Signal-to-noise  ratio,  SNR,  is  the  ratio  of  total  average  signal 
power  to  average  noise  power  as  given  by 

<l}>  (2-53) 


SNR* 


<1  J  *  <<?> 


<*?> 


^  l  n  *  iho> 

Averaging  signal  power  over  time  gives 


<  i 


(2-54) 


where 

<lj>  •  average  signal  power  (A2) 

P,  «■  optical  power  from  aerosol  backscatter  (W) 

Assuming  shot  noise  dominates  other  noise  terms,  then  signal  - to-noise 
ratio  reduces  to 

£\n.’l«h„PJ!(*)dKTJ(tf )  (2-55) 


SNR  * 


2  hBR* 


after  substituting  for  <r,>  and  and  simplifying.  The  term  fV, ,'(/?) 


is  the  aerosol  backscatter  coefficient  at  range  /?.  This  equation  is 
the  ideal  signal- to-noise  ratio  for  an  unfocused  heterodyne  lidar 
assuming  each  scattering  volume  acts  as  a  point  target. 

A  number  researchers  have  studied  receiving  efficiency  for  optical 
heterodyne  receivers.  Sonnenschein  and  Horrigan  (1971)  derived  a  signal 
reduction  factor  to  take  into  account  the  effect  of  wavefront  curvature 
of  the  back9cattered  light  and  the  finite  size  of  the  scattering  volume. 


The  authors  assumed  an  untruncated  Gaussian  output  beam,  spherical  scat- 
terers,  and  a  Gaussian  plane  wave  output  from  the  local  oscillator. 

Clifford  and  Wandzura  (1981)  modified  the  signal  reduction  factor  to 
include  the  reduction  of  spatial  coherence  of  the  backscattered  light 
due  to  atmospheric  refractive  turbulence.  The  resulting  signal  reduc¬ 
tion  factor,  which  multiplies  the  signal- to-noise  ratio  above,  is  given 
by 


SRF  ■ 


(2-56) 


where 


SRF  -  signal  reduction  factor 
D  -  diameter  of  receiver  (m) 
ra  -  transverse  coherence  length  (m) 
f  -  telescope  focus  range  (m) 

Incorporating  the  signal  reduction  factor  and  an  additional  factor, 

.-■I  i,  to  account  for  the  fact  that  the  system  may  not  be  completely  shot 

noise  limited  (Post,  1985),  the  final  form  of  the  heterodyne  lidar 

signal- to-noise  ratio  becomes 

ff.n.n*n,.fl:(#)ni>i*T*(#)/i1s#/?  (2-57) 

SNR  ■=  - - - ; - 

QhBR 7 


where  the  receiving  aperture  is  assumed  to  be  circular  with  area  ii/J'/l, 


2-JL4-. _ Yeloslty..  Ea  llmafc  1  an 

Since  our  primary  objective  is  to  measure  airflow  characteristics 
within  an  aircraft  wake,  we  need  a  method  of  extracting  this  information 
from  the  backscattered  signal.  The  two  most  meaningful  parameters  to 
describe  the  airflow  characteristics  in  a  scattering  volume  are  mean 
velocity  and  velocity  width,  Mean  velocity  is  simply  the  average  veloc¬ 
ity  of  aerosols  along  the  beam  path  within  a  range  gate,  Velocity  width 
indicates  the  variance  or  spread  of  aerosol  velocities  within  a 
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scattering  volume.  Most  of  the  existing  velocity  and  velocity  width 
estimation  techniques  have  been  developed  for  Doppler  weather  radar  and 
are  easily  adapted  for  Doppler  Ildar. 

The  pulse  pair  estimation  algorithm  is  probably  the  most  widely 
used  estimator  of  velocity  and  velocity  width,  Zrnic  (1977)  described 
the  pulse  pair  algorithm  for  estimating  first  three  spectral  moments 
(signal  power,  mean  velocity,  and  velocity  width)  from  Doppler  radar 
weather  echoes.  In  a  later  paper  (Zrnic,  1979)  he  compares  the  pulse 
pair  estimator  to  estimates  derived  from  spectral  processing  using  Four¬ 
ier  transforms  and  then  compares  the  estimation  errors  to  the  Cramer-Rao 
lower  bound  error  derived  from  Maximum  Likelihood  estimates.  It  is  the 
estimation  error  expressions  that  we  are  most  intersted  in  for  our  per¬ 
formance  analysis. 

To  simplify  the  analysis  procedure,  we  will  only  be  interested  in 
estimation  of  velocity  (or  spectral)  width.  The  primary  justification 
for  this  decision  is  that  typical  range  gate  lengths  ate  longer  than  the 
transverse  wake  dimensions  and  thus  the  average  velocity  in  the  range 
gate  generally  averages  to  zero.  In  addition,  we  will  only  look  at  the 
Cramer-Rao  lower  bound  error,  for  width  estimation,  keeping  in  mind  chat 
most  practical  signal  processors  wiit  not  be  able  to  equal  this  perform¬ 
ance.  This  decision  was  made  because  the  error  expressions  for  the 
pulse  pair  and  spectral  processing  width  estimators  are  very  complex  and 
not  really  worth  the  extra  analysis  effort. 

Since  the  Doppler  shift  is  much  higher  at  infrared  wavelengths  chan 
at  microwave  wavelengths,  complete  velocity  information  may  be  extracted 
from  a  single  Ildar  pulse,  as  compared  to  many  radar  pulses.  Thus  the 
Cramer-Rao  error  expression  given  by  Zrnic  was  modified  for  lidar  by 
making  the  following  substitutions  and  redefinitions: 


I 
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T ,  -  T  ~  intrapulse  sample  time 

p  -  aT ,  -  ratio  of  signal  spectral  width  ( w )  to  the  full  electronic 
bandwidth 

tp  -  MT,  -  pulse  duration 
M  -  number  of  intrapulse  samplaa 

In  addition,  tha  axpraasion  has  baan  converted  from  frequency  domain,  as 

given  by  Zrnic,  to  velocity  domain  using 

X  (2-58) 

"  2°" 


whore 

Opi  -  standard  deviation  of  velocity  width  estimate  (m/s) 

Oj,  -  standard  deviation  of  frequency  width  estimate  (Hz) 

So  finally,  the  Cramer-Rao  lower  bound  velocity  width  error  is 
given  by  (Zrnic,  1979:eq.'s  A. 39  and  A, 41) 

for  high  SNR:  (2*59) 

o 

and  for  low  SNR: 


t  „Tt 


These  expressions  indicate  the  typical  errors  associated  with  estimating 
velocity  spread  or  width. 
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Ac  chit  point,  some  discussion  of  preliminary  considerations  is  in 
order  co  sec  the  groundwork  for  analysis  of  atmospheric  and  exhaust 
aerosols.  This  chapter  lays  out  the  criteria  that  will  be  used  in  the 
analysis  along  with  applicable  assumptions  and  justifications.  In  par¬ 
ticular,  lidar  system  parameters  to  be  used  will  be  given,  wake  parame¬ 
ters  for  four  different  test  case  aircraft  will  be  computed,  detection 
geometry  will  be  defined,  signal  processing  requirements  will  be 
established,  and  molecular  absorption  effects  will  be  quantified. 

I~L — Lidar -System  Parameters 

Most  of  the  base  lidAr  system  parameters  to  be  used  in  this  analy¬ 
sis  were  taken  from  Thomson  at  al.  (1989).  Table  3-1  lists  the  base 
system  parameters  for  all  three  Ildars.  The  laser  parameters,  pulse 
energy  and  pulse  width,  reflect  modest  requirements  for  an  airborne 
lidar  system.  These  parameters  do  not  necessarily  imply  the  existence 
of  such  lasers;  however,  active  development  work  on  all  three  lasers  is 
currently  underway,  and  it  is  not  unreasonable  to  assume  that  they  will 
be  available  in  the  next  few  years  (Thomson  ot  al. ,  1989:36,60-62). 

The  detector  quantum  efficiencies  used  were  verified  by  knowledge 
able  researchers.  In  the  9-10  jim  band,  Mercury-Cadmium-Telluride 
(HgCdTe)  detectors  are  available,  both  photoconductive  and  photovoltaic, 
with  quantum  efficiencies  from  60-90%  at  bandwidths  below  300  MHs 
(Schreiber,  1990).  A  nominal  value  of  80%  was  chosen.  For  the  1-2  /urn 
band,  Indium-Gallium-Arsenide  (InGaAs)  and  Germanium  (Ge)  detectors  are 
available  with  quantum  efficiencies  of  70-80%  at  low  bandwidths  (Spears, 
1990).  Again,  a  value  of  80%  was  chosen  for  both  1.064  and  2.091  m 
detectors.  Electronic  bandwidth  for  each  wavelength  corresponds  to  a 
velocity  bandwidth  of  50  m/s . 
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Table  3-1.  Base  Lid&r  System  Parameters 


Parameter 

E939E9 

ES9B9 

Laser  Medium 

Nd :  YAG 

Ho: YAG 

Wavelength,  X  (*«m) 

1.064 

2.091 

9.115 

Pulse  Energy,  B  (J) 

0,2 

0.2 

0.2 

Pulse  Duration,  x„  (^s) 

0.1 

■om 

1.0 

Receive/Transmit  Aperature  Diame¬ 
ter,  D  (m) 

mm 

0.3 

mm 

Optics  Transmission  (two-way),  n. 

0.5 

0.5 

0.5 

Heterodyne  Mixing  Effloieney,  n* 

0.3 

0.5 

Detector  Quantum  Efficiency, 

0.8 

0.8 

0.8 

Shot  Noise  Correction  Term,  A, 

0.5 

0.5 

0.5 

Electronic  Bandwidth,  B  (MHz) 

94 

48 

11 

Telescope  Focus  Range,  /  (km) 

m 

m 

Detection  out  to  30  miles  (80  km)  in  range  is  desirable.  Lldar 
performance,  therefore,  will  be  evaluated  at  ranges  between  0-80  km.  An 
unfocused  telescope  (focused  at  •)  will  be  used  in  all  calculations. 
While  the  signal-to-noise  ratio  at  a  given  range  may  be  Improved 
significantly  by  focusing  the  telescope  at  or  somewhat  farther  than  that 
range,  the  signal-to-noise  ratio  at  other  ranges  suffers  (Thomson  et 
al.,  1989).  In  addition,  for  Gaussian  output  beams,  focusing  at  ranges 
beyond  the  Rayleigh  range  of  the  telescope  is  not  possible  due  to 
diffraction  limitations,  The  Rayleigh  range  is  given  by  the  equation 
(Verdeyen,  1989 :  G9) 

mi'l  (3-1) 

z.-T- 

where 


zQ  -  Rayleigh  range  (m) 

we  -  radius  of  beam  waist  (m) 

\  -  wavelength  (m) 
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The  Rayleigh  ranges  for  the  three  wavelengths  of  interest  and  for  a  0,3 
m  telescope  ere  66.4  kit,  33.8  km,  end  7.7  km  et  1.064  nm,  2.091  Atm,  end 
9.115  A*m,  respectively.  Furthermore,  prectlcel  hardware  considerations 
limit  telescope  focusing  to  less  then  about  30  km  (Zhao  end  Post, 
1989:188). 

1.2 _ _ Aircraft  Examples 

Four  hypothetical  aircraft  will  be  used  to  analyze  wake  detection. 
These  aircraft  were  chosen  to  span  a  wide  range  of  potential  aircraft 
sices  and  missions.  Included  are  large  and  medium  transport  aircraft 
and  large  and  small  fighter  aircraft.  The  pertinent  dimensions  and 
cruise  parameters  for  each  of  these  aircraft  are  listed  in  Table  3*2. 
Note  that  aircraft  velocities  used  here  are  subsonic  or  transonic. 
Supersonic  velocities  will  not  be  considered. 

Table  3-2.  Test  Aircraft  Dimensions  and  Cruiss  Parameters 


Parameter 

Large 

Transport 

Medium 

Transport 

Large 

Fighter 

Small 

Fighter 

Wingspan,  s  (m) 

63 

45 

13 

8 

Mass  (fully  loaded),  Af 
(kg) 

350,000 

135,000 

25,000 

11,000 

Cruise  Speed,  1'0  (m/s) 

250 

235 

300 

300 

Cruise  Altitude,  h  (km) 

10.3 

10.3 

12.2 

12.2 

Air  Density  at  Altitude, 

P  (kg/ms) 

0.398 

0.398 

0.302 

0.302 

Wing  Loading  Constant,  K 

0,785 

0.783 

0.785 

0.785 

Engines 

4  (a  40,000 
lbs  thrust 
each 

4  @  14,000 
lbs  thrust 

each 

2  <§ 

25,000 

lbs 

thrust 

each 

1  @ 

10,000 

lbs 

thrust 

Applying  the  wake  models  presented  in  chapter  2,  we  may  calculate 
the  size  of  the  stable  wake  and  the  air  velocities  within  it  for  each  of 
these  aircraft  at  the  given  cruise  conditions.  Table  3-3  lists  these 
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Tangential  Velocity  (m/s) 


Figure  3-1.  Tangential  Velocity  Profiles  of  Test  Aircraft 


wake  parameters  for  all  four  aircraft.  In  addition,  Figure  3-1  shows 
the  tangential  air  velocity  about  one  of  the  wake  vortices  as  a  function 
of  distance  from  the  vortex  center. 

Table  3-3.  Wake  Parameters  for  the  Four  Test  Aircraft 


Parameter 


Vortex  Separation,  6  (m) 


Vortex  Core  Radius,  re 
(m) 


Peak  Tangential  Velocity, 

k  /,*n k  (it/ ®  ) 


Sink  Velocity,  I',  (m/s) 


Core  Drag  Velocity,  I’, 
(m/a) 


Thrust  Velocity,  V,  (m/s) 


Recirculation  Cell  Width 
(m) 


Recirculation  Call  Haight 
<») 


Wake  Length,  /.„  (km) 


To  reduce  the  analysis  problem  to  a  manageable  level,  we  must  make 
some  assumptions  about  the  detection  geomotry.  First,  let  us  assume 
that  the  laser  beam  passes  horizontally  through  the  wake.  This  assump¬ 
tion  is  a  good  one  since  we  are  primarily  interested  in  near  horizontal 
air-to-air  detection  at  long  range.  Second,  we  may  assume  that  the 
transverse  dimension  of  the  laser  beam  (the  width)  is  small  compared  to 
the  wake  dimensions.  Assuming  the  transmitted  beam  has  a  Gaussian  pro¬ 
file  and  is  unfocused,  the  beam  spot  size  at  a  distance  n  from  the  lidar 
transmitter  is  given  by  (Verdeyen,  1989:69) 

T  l  9  \,T/*  0*2) 


z )  •  w 


•['*(r.)7 
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where 


tu(z)  -  beam  apoc  radius  at  distance  z  (m) 

2,  -  Rayleigh  range  (m) 

Since  the  output  beam  is  roughly  0.3  m  in  diameter  (see  Table  3*1), 
the  beam  diameters  at  a  range  of  80  km  are  0.47  m,  0.77  m,  and  3.1  m  at 
1.064  pm,  2.091  pm,  and  9.115  /im,  respectively.  These  beams  are  all 
smaller  than  the  recirculation  call  diameters  given  in  Table  3-3. 

In  analyzing  return  signals  from  aircraft  wakes,  we  therefore  only 
need  to  consider  the  vertical  height  of  the  beam  within  the  wake  and  the 
horizontal  angle  between  the  laser  and  the  aircraft  flight  direction. 
Figure  3*2  shows  this  simplified  geometry. 

Two  final  geometric  assumptions  concerning  the  beam  path  will  be 
made.  First,  we  will  use  the  flat  earth  assumption;  that  is,  curvature 
of  the  earth  will  be  neglected.  This  assumption  leads  to  an  error  in 
vertical  height  of  about  0.2  km  over  a  50  km  horizontal  path  on  the 
ground.  Second,  we  will  neglect  any  refractive  bending  of  the  lidnr 
beam  due  to  the  change  in  atmospheric  index  of  refraction  with  altitude. 
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Target  Aircraft 


Figure  3*2.  Simplified  Wake  Detection  Geometry 
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The  return  signal  from  any  given  range  gate  interval  contains  spec¬ 
tral  information  on  the  distribution  of  aerosol  velocities  within  the 
scattering  volume.  In  order  to  better  understand  how  to  extract  this 
information  and  Identify  signals  containing  wake  signatures,  we  should 
see  what  the  returns  from  a  wake  look  like. 

To  get  an  idea  of  the  general  shape  of  the  return  spectrum  from 
within  a  wake,  a  computer  simulation  was  performed  to  generate  a  two 
dimensional  flow  field  in  the  cross  section  of  the  medium  transport 
wake.  The  return  signal  consisted  of  the  sum  of  all  the  horizontal 
velocity  components  along  a  150  m  segment,  or  "range  gate",  cantered 
horizontally  on  the  wake.  The  "beam"  was  then  scanned  vertically  from 
the  vortex  center  to  see  how  the  return  spectrum  varied  with  vertical 
dlstanoe.  The  simulation  assumsd  that  the  beam  width  was  smaller  than 
the  2  m  resolution  of  sample  points  within  the  flow  field.  In  addition 
the  lidar  beam  was  assumed  to  be  perpendicular  to  the  length  of  the 
wake.  Figure  3-3  show  these  simulated  return  spectrums  for  vertical 
samples  from  2-14  m  from  the  vortex  center. 

As  Figure  3-3  indicates,  the  return  signal  from  the  wake  center  is 
very  narrow,  with  all  the  energy  centered  about  the  zero  velocity  point 
As  thn  beam  moves  vertically  up  or  down  from  the  center,  however,  the 
return  spectrum  quickly  becomes  very  broad,  with  significant  energy  in 
the  higher  velocity  components,  The  majority  of  the  signal  energy 
remains,  though,  near  zero  velocity.  As  the  beam  continues  to  scan 
vertically  up  or  down,  th#  spectrum  slowly  begins  to  narrow  again, 

While  the  spectra  do  only  have  one  main  lobe,  they  do  not  appear  to  be 
Gaussian  in  shaps.  In  all  cases,  since  the  simulated  range  gate  was 
centered  on  the  wake,  the  average  velocity  was  zero. 
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Figure 


3*3,  Simulated  Return  Spectra  for  Several  Vertical  Snmpl 
through  tha  Wake  of  a  Medium  Sized  Transport  Aircraft 


Velocity  Width  (m/s) 


Vertical  Distance  from  Vortex  Center  (m) 


Figure  3-4,  Velocity  Width  of  Return  Spectra  from  the  Vertically 
Scanned  Wakes  of  U  Aircraft 


3-10 


Probably  the  most  obvious  feature  of  the  return  spectra  from  the 

wake  Is  the  wide  spread  of  velocities.  Figure  3-4  shows  the  velocity 
width  (taken  as  twice  the  velocity  standard  deviation)  of  the  return 
signal  plotted  versus  vertical  distance  from  the  vortex  center  for  the 
four  test  aircraft.  The  graph  clearly  shows  the  wide  spread  in  velocity 
distribution  near  the  vortex  cores.  Since  typical  spectrum  width  due  to 
atmospheric  turbidity  is  only  about  1-2  m/s,  the  6-12  m/s  spread  due  to 
wake  turbulence  seems  to  be  a  good  indicator  of  the  presence  of  a  wake 
in  the  return  signal. 

Range  gate  width  affectR  the  return  spectrum  a  great  deal.  For  the 
previous  simulation,  the  range  gate  interval  was  150  m,  which  corre¬ 
sponds  to  a  pulse  duration  of  1  /js.  Better  range  resolution  may  be 
obtained  by  using  shorter  pulse  durations,  but  this  improvement  comes  at 
the  expense  of  lowe'‘  frequency  resolution.  On  the  other  hand,  better 
frequency/velocity  - . solution  may  be  obtained  at  the  expense  of  longer 
range  resolution  by  using  longer  pulses.  Figures  3-5. a-d  show  the 
effect  of  different  range  gate  lengths  (80,  150,  and  750  m)  on  the 
return  signal  velocity  width  for  the  four  test  aircraft.  As  the  graphs 
show,  the  velocity  width  decreases  with  increasing  range  gate  length 
since  the  turbulent  ake  makes  up  a  smaller  fraction  of  the  total  sig¬ 
nal  . 

To  estimate  how  much  s ignal - to-noise  is  required  to  detect  w;  ke 
signals  we  will  use  the  Cramer-Rno  lower  bound  error  for  the  Maximum 
Likelihood  velocity  width  estimator  as  presented  in  chapter  2.  Plots  of 
width  estimation  error  are  shown  in  Figures  3-6. a  and  3-6. b  for  four 
different  "true"  spectral  widths.  Figure  3-6. a  is  for  the  1.064  case 
with  a  pulse  duration  of  0.1  ju s,  while  Figure  3-6.b  is  for  the  same 
wavelength  with  a  0.5  jis  pulse  duration.  Due  to  wavelength  scaling 
effects,  these  plot3  are  approximately  identical  to  tho  error  plots  for 
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Figure  3*5. a.  Velocity  Width  of  Return  Spectra  from  the  Large  Transport 
Wake  for  3  Different  Range  Gate  Intervals 
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5,b,  Velocity  Width  of  Return  Spectra  from  the  Medium  Trans¬ 
port  Wake  for  3  Different  Range  Gate  Intervals 
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Velocity  Width  (m/s) 


Figure  3-5.c.  Velocity  Width  of  Return  Spectra  from  the  Large  Fighter 
Wake  for  3  Different  Range  Gate  Intervals 
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Figure  3-5, d.  Velocity  Width  of  Return  Spectra  from  the  Small  Fighter 
Wake  for  3  Different  Range  Gate  Intervals 
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Signal— to— -Noise  Ratio 


Figure  3*6.a,  Velocity  Width  Estimation  Error  for  Short  Pulse  Duration 

and  Four  "True"  Spectral  Widths 
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Figure  3-6 


Signal-to-Noise  Ratio 


.b.  Velocity  Width  Estimation  Error  for  Long  Pulse  Duration 
ana  Four  "True"  Spectral  Widths 
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0.2  /js  and  1.0  /»s  pulse  durations  at  2.091  and  1.0  /is  and  5.0  /is 
pulse  durations  at  9.115  /4m.  The  plots  were  generated  assuming  Nyquist 
sampling  and  a  velocity  bandwidth  of  50  m/a. 

Figure  3*6  shows  that  estimation  error  Increases  with  increasing 
spectral  width.  In  addition,  estimation  error  is  higher  for  the  shorter 
pulse  length  case  since  the  shorter  range  gate  interval  contains  fewer 
signal  samples.  Minimum  "reasonable"  signal- to-noise  ratios  are  in  the 
vicinity  of  0  to  +5  dB  for  the  short  pulse  length  case  and  -5  to  0  dB 
for  the  long  pulse  length  case.  Since  the  Cramer-Rao  error  is  u  lower 
bound,  we  will  use  the  more  pessimistic  values  of  3  dB  and  0  dB  for  the 
short  and  long  pulse  length  cases,  respectively.  At  these  signal-to- 
noise  ratios,  the  optimum  width  estimators  are  within  50%  of  the  true 
width  for  velocity  widths  below  15  m/s.  Thus,  signals  below  0  dB  in 
wideband  power  signal- to-noise  ratio  are  not  very  useful  for  wake  detec¬ 
tion  using  the  spectral  width  method. 

_ Molecular  Absorption  Effects 

Three  runs  of  the  computer  program  FASCODE  were  made  to  estimate 
the  magnitude  of  molecular  absorption  in  the  atmosphere  at  the  three 
laser  wavelengths.  The  geometry  used  was  chosen  to  give  a  worst  case 
scenario  to  see  if  molecular  absorption  was  a  factor  or  not.  The  pro¬ 
gram  computed  the  one-way  transmission  over  a  50  kin  long  path  from  an 
initial  height  of  12.2  km  (40,000  ft)  to  a  final  height  of  7.62  km 
(25,000  ft).  The  resulting  look  angle  was  about  5  degrees  down.  Table 
3-4  lists  the  resulting  one-  and  two-way  transmissions  As  Table  3-4 
clearly  ohows ,  molecular  absorption  along  the  beam  path  is  negligible 
for  this  analysis. 
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Table  3-4 


Atmospheric  Transmission  due  to  Molecular  Absorp 
50  km  Path  from  12.2  km  to  7.62  km  Altitude 


Transmission 

Wavelength 

One • way 

Two - way 

1,064  pm 

0.9994 

0.9987 

2.091  pm 

0.9906 

0.9813 

9.115  Min 

0.9934 

0.9868 

ion  for  a 
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. A  t  rnos  phexic..  Ac  ros  o  la 


In  this  chnpter,  wo  will  try  to  ooe  what  effect,  atmospheric  aero¬ 
sols  have  on  the  performance  of  a  Doppler  litlnr  wake  detection  system. 
Atmospheric  aerosols,  at  least  in  this  report,  refer  to  naturally 
occurring  aerosol  particles.  Dense  aerosol  regions ,  such  as  clouds, 
will  not  he  analyzed  since  it  is  assumed  that  Ildar  beams  cannot  probe 
very  far  into  them.  We  will  assume  a  cloud-free  line  of  sight  (CFLOS). 

This  chapter  is  broken  down  into  throe  main  subsections.  In  the 
flrnt  subsection,  the  vortical  structure  of  aerosols  in  the  atmosphere 
will  bo  described,  In  the  second  subsection,  models  of  aerosol  back- 
scatter  and  extinction  coefficients  will  be  presented  for  use  an  test 
cases.  Finally,  in  the  third  subsection,  signal- to-noise  ratio  will  be 
analyzed  for  those  model  atmospheres.  Wake  detectability  assuming  scat¬ 
tering  from  atmospheric  aerosols  alone  will  bo  assessed. 

ul  .AtraxoL  JLtxuc  £  uru 

The  atmosphere  consists  of  several  well  defined  altitude  regions, 
each  with  its  own  unique  properties,  below  30  km  in  altitude,  these 
regions  Include  tho  boundary  layer,  the  free  troposphere ,  and  the 
stratosphere.  Idaeh  region  is  described  below. 

(t . .l.i  1  .-..-JUyumltu'y.. Layer 

The  boundary  layer  is  the  lower  part  ol  the  troposphere  and  exists 
from  tho  earth's  suriaco  up  to  around  2  km  In  height.  The  salient 
feature  of  the  boundary  layer  is  u  very  high  aerosol  content  (the  haze 
wc  see  as  wo  look  toward  the  horizon).  In  addition,  the  boundary  layer 
aerosols  tend  to  be  larger  than  aerosols  at  higher  altitudes  (heavier 
particles  sinks).  Thus,  bnckscatter  anti  attenuation  in  the  boundary 
layer  are  generally  very  high.  Boundary  layer  aerosol  typo  depends 
largely  on  geographic  location. 
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For  this  study,  boundary  layer  aerosols  are  not  very  important 
since  we  ere  only  interested  in  detecting  aircraft  at  cruise  altitudes. 

4.1.2.  Free  Troposphere 

The  free  troposphere  is  the  upper  portion  of  the  troposphere  end 
extends  from  the  top  of  the  boundary  layer  to  the  tropopause,  which 
separates  the  troposphere  from  the  statosphere.  Tropopause  height 
varies  both  seasonally  and  geographically,  but  is  generally  between  8*17 
km,  Mean  tropopause  height  is  highest  at  latitudes  near  the  equator  and 
lowest  at  latitudes  near  the  poles.  Daily  variations  in  tropopause 
height  often  exceed  seasonal  variations  (Brown  and  Kunkel,  1985:16-46). 
Aerosol  concentrations  in  the  free  troposphere  decrease  rapidly  from  the 
boundary  layer  up  to  the  tropopause.  In  contrast  to  the  boundary  layer, 
the  tropopause  is  cleanest  part  of  the  lower  atmosphere  and  thus  exhib¬ 
its  relatively  low  backacatter.  Aerosols  in  the  free  troposphere  gener¬ 
ally  consist  of  water  soluble  particles  (ammonium  and  calcium  sulfate 
and  organic  compounds)  and  dust-like  particles  (Fenn  et  al,, 

19B5 : 18-12) . 

Aerosol  concentration  and  type  in  the  free  troposphere  are  of  great 
Importance  in  this  study  since  many  aircraft  cruise  at  these  altitudes. 
Since  this  is  typically  the  most  aerosol-free  layer  in  the  lower  atmo¬ 
sphere,  it  is  expected  that  wake  detection  here  will  be  the  most  diffi¬ 
cult. 

ft  .1.3. . -S.,t£At.Q.apharj8 

Above  the  troposphere  lies  the  stratosphere.  The  stratosphere 
extends  from  the  top  of  the  tropopause  up  to  around  45  km  in  altitude. 
For  this  analysis,  we  are  only  interested  in  altitudes  up  to  about  20  km 
(roughly  60,000  ft),  which  is  the  limit  of  most  aircraft.  The  majority 
of  aerosols  in  the  stratosphere  are  of  volcanic  origin.  During  times  of 
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Increased  volcanic  activity,  the  stratosphere  becomes  heavily  laden  with 
volcanic  debris  and  thus  exhibits  very  high  backscatter,  Over  time, 
these  aerosols  tend  to  fall  out,  and  the  stratosphere  becomes  slowly 
cleaner.  There  is  some  evidence  that  suggests  that  volcanic  aerosols 
tend  to  filter  out  only  at  the  higher  latitudes  near  the  earth's  poles 
since  tropical  and  mid* latitude  tropopauses  appear  to  remain  very  clean 
despite  high  volcanic  loading  in  the  stratosphere  (Kent  et  al., 
1986:2-31). 

Most  military  aircraft  ars  capable  of  flying  at  stratospheric  alti¬ 
tudes,  For  this  reason,  the  stratosphere  is  an  important  region  in  this 
study.  Since  stratospheric  aerosol  loading  can  vary  a  great  deal 
depending  on  recent  volcanic  activity,  it  is  expected  that  wake  detect¬ 
ability  in  the  stratosphere  will  also  vary  greatly. 

LJL. _ Backscatter-  and  Extinction 

The  general  atmospheric  structure  described  above  is  fairly  well 
known,  The  problem  is  in  quantifying  what  the  backscatter  and  extinc¬ 
tion  levels  are  for  a  given  wavelength  of  light,  at  a  given  altitude,  at 
a  given  geographic  location,  and  during  a  given  season.  Furthermore,  we 
would  like  to  know  how  much  variability  exists  in  thess  levels  and  how 
often  they  occur,  Unfortunately,  a  complete  statistical  modal  of  global 
backscatter  and  oxtinctlon  does  not  exist.  Adequate  models  of  mean 
aerosol  characteristics  are  available  end  they  will  be  described  below. 

Using  these  aerosol  models,  three  test  atmospheres  were  developed 
to  evaluate  Ildar  performance  under  a  wide  range  of  atmospheric  aerosol 
conditions.  In  this  analysis,  the  atmosphere  is  assumed  to  be  horizon¬ 
tally  homogeneous;  that  is,  backscatter  and  extinction  are  constant  at 
any  given  altitude, 
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Finally,  absorption  of  laser  light  by  atmospheric  gas  molecules 
must  not  be  overlooked.  The  results  of  computer  analyses  for  the  three 
wavelengths  of  Interest  are  presented. 

4-2.1. _ Geophysics  Laboratory  Aerosol  Models 

The  Air  Force's  Geophysics  Laboratory  has  developed  an  extensive 
set  of  atmospheric  aerosol  models  (Fenn  et  al.,  1985: section  18-2).  The 
models  consist  of  particle  size  distributions  for  several  types  of  aero* 
sol  in  the  boundary  layer,  free  troposphere,  and  stratosphere,  The  size 
distributions  for  each  aerosol  type  are  baeed  on  experimental 
measurements  by  many  researchers  and  assume  spherical  particles.  Opti¬ 
cal  scattering  and  absorption  parameters  are  calculated  from  these  mod¬ 
els  using  Mie  scattering  theory.  Variations  in  the  aerosol 
characteristics  due  to  humidity,  visibility,  geographic  location, 
season,  and  volcanic  aotivlty  are  included  in  the  models. 

The  Geophysics  Laboratory  aerosol  models  have  been  Implemented  on  a 
computer  program  called  BACKSCAT  that  simulates  lidar  backseatter  for  a 
wide  range  of  wavelengths  from  0.20  to  40.0  pm  (Gutvons  et  al.,  1988). 
BACKSCAT  runs  on  an  IBM- PC  type  computer  and  computes  backscattered 
power  versus  range  for  any  given  Ildar  system  and  atmospheric  model.  In 
addition,  the  program  generates  altitude  profiles  of  aerosol  and  molecu¬ 
lar  backseatter  and  extinction. 

While  these  models  do  give  mean  backseatter  and  extinction  levels 
for  different  atmoepheres,  they  do  not  indicate  how  much  variability  may 
occur  in  the  levels.  Furthermore,  the  models  are  based  on  assumed  size 
distributions  end  aerosol  types  and  not  on  actual  lidar  measurements. 

As  discussed  later,  some  limited  validation  of  the  models,  st  least  at 
0.5  pm  and  10,6  pm,  indicates  good  correlation  with  lidar  data  (Koenig, 
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1990).  Though  not  ideal,  the  Geophysics  Laboratory  aerosol  models 
appear  to  be  one  of  the  most  extensive  available  and  so  will  be  used  in 

this  study  to  generate  model  atmospheres  to  use  as  test  cases. 

4.2.2.  Modal  Atmospheres  for  Calculations 

BACKSCAT  was  used  to  generate  three  model  atmospheres  for  use  as 
test  cases  to  compute  signal-to-noisa  ratio.  The  three  teat  cases  rep¬ 
resent  atmospheres  with  a  wide  range  of  aerosol  loading  from  very  clean 
to  very  dirty  air.  The  atmospheric  modols  consist  of  altitude  profiles 
of  backscatter  and  extinction  coefficients  from  both  aerosols  and  mole¬ 
cules  at  the  three  infrared  wavelengths  of  interest.  In  all  cases,  the 
tropopause  height  was  9  km,  the  top  of  the  stratosphere  was  at  29  kin, 
and  the  top  of  the  boundary  layer  was  at  2  km.  These  values  were  fixed 
by  BACKSCAT.  Hopefully,  future  upgrades  to  this  program  will  allow 
layer  heights  to  be  changed. 

The  molecular  backscatter  and  total  scattering  coefficients  were 
computed  by  BACKSCAT  using  Che  Rayleigh  scattering  theory  presented  in 
ohapter  2,  These  parameters  do  not  change  with  varying  aerosol  models. 
Molecular  backscatter  coefficients  for  the  three  wavelengths  are  shown 
in  Figure  4-1. a  and  total  scattering  coefficients  are  shown  in  Figure 
4-l.b.  For  negligibly  small  molecular  absorption,  the  total  molecular 
scattering  coefficients  are  also  the  molecular  extinction  coefficients. 

The  clean  atmosphere  was  developed  to  give  the  lowest  possible 
backscatter  and  extinction  values  at  all  altitudos.  The  boundary  layer 
model  used  wee  the  rural  aerosol  model  at  01  humidity  and  50  km  visibil¬ 
ity.  The  tropospheric  model  also  used  01  relative  humidity.  For  the 
stratosphere ,  the  background  fall/winter  model  was  Implemented,  which 
assumed  almost  no  volcanic  aerosol  loading.  Figure  4- 2. a  shows  the 
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Figure  4-2. b.  Aerosol  Extinction  Coefficient  Profile  for  the  Clean 

Atmosphere  Model 
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Figure  4*3. e.  Aerosol  Backscatter  Coefficient  Profile  for  the  Moderate 

Atmosphere  Model 
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Figure  4,3.b.  Aerosol  Extinction  Coefficient  Profile  for  the  Moderate 

Atmosphere  Model 


resulting  aerosol  backscatter  coefficient  profiles  and  Figure  4*2. b 
shows  the  aerosol  extinction  coefficient  profiles  for  the  clean  atmo¬ 
sphere. 

The  moderate  atmosphere  was  developed  to  be  give  backscatter  values 
roughly  halfway  between  the  clean  and  dirty  atmospheres.  In  no  way, 
however,  does  this  imply  that  the  moderate  atmosphere  is  representative 
of  statistical  average  or  moat  probable  backscatter  levels.  The  bound¬ 
ary  layer  model  used  was  a  rural  aerosol  model  at  70%  relative  humidity 
and  23  km  visibility.  The  tropospheric  humidity  was  also  set  at  70%, 
and  the  stratospheric  model  used  was  the  spring/summer  aged-moderate 
volcanic  model.  The  aged  volcanic  model  assumes  that  enough  time  has 
elapsed  since  injection  that  the  larger  volcanic  particles  have  settled 
out.  The  moderate  volcanic  model  relates  to  the  total  particle  number 
density.  Figure  4- 3. a  shows  the  resulting  aerosol  backscatter  coeffi¬ 
cient  profiles  for  tho  moderate  atmosphere.  Figure  4-3. b  shows  the 
resulting  aerosol  extinction  coefficient  profiles. 

The  dirty  atmosphere  model  was  dsveloped  to  give  the  highest  possi¬ 
ble  backscatter  values  ovsr  all  altitudes.  The  boundary  layar  model 
used  for  the  dirty  atmosphere  was  the  urban  aerosol  model  at  99% 
relative  humidity  and  2  km  visibility.  In  the  troposphere,  relative 
humidity  was  also  set  to  99%.  The  stratospheric  model  used  the 
spring/summer  fresh-extreme  volcanic  aerosol  model.  The  fresh  volcanic 
model  assumes  a  recent  eruption,  resulting  in  a  lot  of  large  particles. 
Figures  4-4. a-b  show  the  aerosol  backscatter  and  extinction  coefficient 
profiles  of  the  dirty  atmosphere  model. 

In  the  three  atmosphere  models  above,  it  is  apparent  that  the  back* 
scatter  coefficients  for  the  shorter  wavelengths  are  higher  than  those 
for  the  longer  wavelengths.  Also  note  that  molecular  backscatter  Is 
very  significant  at  1.064  pm  compared  to  aerosol  backscatter ,  but  it  is 
completely  negligible  at  9.115  pm.  From  clean  atmosphere  to  dirty 
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atmosphere,  backscatter  in  the  free  troposphere  changes  very  little, 
while  backscatter  in  the  stratosphere  changes  3*4  orders  of  magnitude l 
Boundary  layer  backscatter  varies  by  about  an  order  of  magnitude  from 
clean  to  dirty. 

4.2.3, _ Comparison  with  Measured  Aaroaol  Data 

The  three  models  of  atmospheric  aerosol  backscatter  and  extinction 
ware  compared  to  measured  data  for  validation,  Unfortunately,  no  usable 
lidar  data  was  found  for  the  1.064  and  2.091  /am  cases;  however,  some 
published  1  urn  extinction  data  from  the  SAGE  1  satellite  was  located  for 
comparison.  Figure  4*5  shows  two  SAGE  I  extinction  coefficient  pro* 
files.  One  profile,  taken  between  September*November  1979,  represents  a 
relatively  clean  atmosphere,  free  of  volcanic  aerosols  (Kent  et  al., 
1966:2*17),  Th«  second  profile,  taken  between  June-August  1980,  shows 
the  effects  of  the  eruption  of  Mount  St.  Helens  (Kent  et  al., 

1986:2-29).  Both  profiles  were  taken  between  20°-40B  north  latitude. 

The  extinction  data  in  Figure  4-5  compares  well  with  the  clean  and 
moderate  atmosphere  models;  however,  the  dirty  atmosphere  model  shows 
much  more  volcanic  enhancement  chan  the  SAGE  I  data  following  the 
eruption  of  Mount  St.  Helens.  For  this  reason,  the  1.064  nm  dirty 
atmosphere  model  should  be  taken  as  an  absolute  extreme  case. 

A  fair  amount  of  published  10,6  jam  lidar  data  exists,  which  was 
used  for  comparison  with  the  9,115  models.  Backscatter  coefficients 
in  the  clean  atmosphere  model  compare  quite  well  with  10.6  #am  lidar  data 
measured  in  the  northern  hemisphere.  Rothermel  et  al.  (1989)  report  "a 
high  frequency  of  occurrence  for  a  narrow  range  of  low  backscatter  val¬ 
ues  in  the  troposphere  •  a  background  mode."  At  altitudes  between  9-12 
km,  the  background  backscatter  coefficient  is  roughly  3xlO*lJ  m"lsr*1 
(Rothermel,  et  al.,  1989:1041).  The  clean  atmosphere  model  developed 
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Figure  4-5.  1  nm  Aerosol  Extinction  Coefficient  Profiles  from  SAGE  1 

Satellite  Data  for  Both  Clean  and  High  Volcanic  Periods 
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above  indicates  backscatter  coefficients  around  5-6xlO*n  m^sr*1  at 
these  altitudes.  The  small  difference  may  be  attributed  to  the  expected 
higher  backscatter  at  9  n  (Ancellet  at  al. ,  1988). 

Volcanic  enhancement  at  9.115  nm  was  also  compared  to  10.6  *xm  lidar 
data.  Post  (1984)  reports  an  increase  in  backscatter  between  about 
18*23  km  altitude  over  Boulder,  CO  following  the  eruption  of  El  Chichon 
in  Mexico  in  1982.  Peak  volcanic  backscatter  coefficient  was  8xl0*9 
m*lsr*1,  which  was  about  an  order  of  magnitude  above  ambient  levels 
before  the  eruption  (Post,  1984:2509).  This  level  of  volcanic  backseat* 
ter  enhancement  is  a  factor  of  8  above  the  1x10*®  nr*ar*1  level  shown  in 
the  moderate  atmosphere  model  but  is  a  almost  two  orders  of  magnitude 
below  the  7xlO*7  m*lar‘l  level  shown  in  the  dirty  atmosphere  model, 

Again,  volcanic  enhancement  in  the  dirty  atmosphere  model  appears  to  be 
at  least  an  order  of  magnitude  too  extreme. 

The  clean  and  moderate  atmosphere  models  compare  well  with  measured 
data  at  1  and  10  /xm,  Stratospheric  backscatter  and  extinction  in  the 
dirty  atmosphere  model  are  1*2  orders  of  magnitude  above  published  val* 
uea.  Results  from  this  model,  therefore,  should  be  accepted  as  absolute 
extremes  at  best. 

LJL — P.ftr£Qrrcance...Aaalyala 

The  results  of  slgnal*to*noise  computations  using  the  model  atmo¬ 
spheres  are  presented  below.  An  analysis  of  these  results  identifies 
implications  for  wake  detectability  from  atmospheric  aerosols  alone. 

UjuL _ Signal  •  so  £annmtailflni 

Signal • to -noise  ratio  was  computed  using  the  equations  developed  in 
chapter  2,  the  lidar  system  parameters  given  in  chapter  3,  and  the  test 
atmospheres  developed  in  this  chapter,  Simulations  were  run  for  two 
different  aircraft  altitudes  (lidar  altitudes):  9  km  (30  kft)  and  13.7 
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km  (45  kfc) .  The  9  km  altitude  is  in  the  tropopause.  Signal- to-noise 
ratio  was  computed  at  each  altitude  for  lidar  beam  angles  from  +6 
degrees  up  to  *6  degrees  down  in  2  degree  increments.  Range  varied  from 
2  to  80  km  in  2  km  increments.  Signel-to-noise  ratio  values  presented 
here  represent  wideband  signal-to-noise  for  a  single  laser  pulse  (no 
pulse  integration) . 

For  each  scenario,  the  signal-to-noise  ratio  versus  range  curves 
for  the  7  look  angles  were  plotted.  The  data  on  these  graphs  was  then 
converted  to  contour  lines  of  constant  signal  level  and  plotted  on 
graphs  of  altitude  versus  horizontal  range.  These  graphs  gave  the 
clearest  picture  of  which  areas  of  the  atmosphere  the  signal  level  was 
sufficient  for  wake  detection,  The  resulting  graphs  are  shown  below  in 
Figures  4-6  through  4-11.  The  cross  hatched  arses  show  where  signal-to- 
noise  ratio  is  above  1  (0  dB).  Contours  are  cubic  spline  fits  to  the 
data  points  shown  on  the  graphs. 

(t  _ Ana  ly  ill  -aUSiiuila 

The  overall  results  for  atmospheric  aerosols  slone  are  rather  dis¬ 
couraging.  In  the  clean  atmosphere  case,  none  of  the  lidars  ahow  any 
appreciable  detection  capability.  Only  the  9.115  pm  Ildar  has  areas 
where  the  signal-to-noise  ratio  exceeded  0  dB,  and  these  are  only  at 
close  range  or  at  very  low  altitudes  where  aerosol  loading  was  high. 
Furthermore,  if  we  use  the  5  dB  minimum  for  short  pulse  lengths  (see 
chapter  3),  the  situation  is  even  worse.  Overall,  wake  detection  for 
the  clean  atmosphere  case  and  the  given  lidar  configurations  is  out  of 
the  question. 

The  results  for  the  moderate  atmosphere  case  are  only  slightly  bet¬ 
ter.  As  before,  the  1.064  pm  and  2.091  pm  lidars  do  not  show  any  real 
wake  detection  capability.  The  9.115  pm  lidar  does,  however,  have  some 
rather  large  areas  where  signal-to-noise  is  over  0  dB.  At  9  km 


4-18 


Clean  Atmoephere 
1.08  pm 
9  km  Altitude 


Horizontal  Range  (km) 


4-6. a.  Signal*to-Noi»e  Ratio  for  Clean  Atmoephere,  Aircraft  at  9 
km  Altitude,  and  1.064 
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Figure  4-6, b.  Signal- to-Noise  Ratio  for  Clean  Atmosphere,  Aircraft  at  9 

km  Altitude,  and  2.091  jjm 
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Figure  4-6. c.  Signal -to-Noise  Ratio  for  Clean  Atmosphere,  Aircraft  at  9 

km  Altitude,  and  9.115  moi 
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Figure  4-7. a.  Signal ■ to -Noise  Ratio  for  Clean  Atmosphere,  Aircraft  at 

13.7  km  Altitude,  and  1.064  *im 
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Figure  4-7, b,  Signal • to -Noise  Ratio  for  Clean  Atmosphere,  Aircraft  at 

13.7  km  Altitude,  and  2.091 
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Figure  4-7. c.  Signal • to-Nolae  Ratio  for  Clean  Atmosphere,  Aircraft  at 

13.7  km  Altitude,  and  9.115  urn 
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Figure  4-8. a.  Signal-to-Noi.se  Ratio  for  Moderate  Atmosphere,  Aircraft 

at  9  km  Altitude,  and  1,064  jim 
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Figure  4*8, b,  Signal* to-Noise  Ratio  for  Moderate  Atmosphere,  Aircrn 

at  9  km  Altitude,  and  2.091  pm 
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Figure  4 • 8 . c .  Signal- to-Noise  Ratio  for  Moderate  Atmosphere, 

at  9  km  Altitude,  and  9,115  urn 
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9. a.  Signal- to-Noist'  Ratio  for  Moderate  Atmosphere, 
at  13.7  km  Altitude,  and  1.064  pm 
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4-9. c,  Signal- to-Noise  Ratio  for  Moderate  Atmosphere, 
at  13.7  km  Altitude,  and  9,115  *im 
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Figure  4-10. a.  Signal - to-Noise  Ratio  for  Dirty  Atmosphere,  Aircraft  at 

9  km  Altitude,  and  1.064  jim 
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Fiture  4-10. c.  Signal-to-Noise  Ratio  for  Dirty  Atmosphere,  Aircraft  at 

9  km  Altitude,  and  9.115  pm 
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tO’Noise  Ratio  for  Dirty  Atmosphere,  Aircraft  at 
.7  km  Altitude,  and  1.064  /urn 
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.b.  Signal-to-Noise  Ratio  for  Dirty  Atmosphere,  Aire 
13.7  km  Altitude,  and  2.091 
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-11. c.  Signal-to-Noise  Ratio  for  Dirty  Atmosphere,  Aircraft  at 
13.7  km  Altitude,  and  9,115  /ini 
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altitude,  maximum  detection  range  is  about  22  km  horizontally,  45  km  at 
+6  degrees  up,  and  over  80  km  at  -6  degrees  down.  As  expected,  the 

troposphere  has  the  lowest  signal-to-noise  values.  At  13.7  km  altitude, 
a  maximum  detection  range  of  about  47  km  is  evident  for  near  horizontal 
look  angles. 

In  the  dirty  atmosphere  case,  the  1.064  pm  Ildar  again  fails  to 
show  any  reasonable  performance.  The  2.091  pm  Ildar  does  exhibit  good 
signal  lavels  from  the  stratospheric  regions  due  to  the  Increased  volca¬ 
nic  aerosol.  Detection  ranges  beyond  80  km  are  evident  for  the  9  km 
altitude  case  looking  up  into  the  volcanic  layer.  At  13.7  km  altitude, 
the  2,091  pm  lidar  shows  maximum  detection  range  from  30  km  looking  down 
to  nearly  80  km  looking  up.  Except  for  the  tropopause  and  free  tropo¬ 
spheric  regions,  the  9.115  pm  lidar  shows  reasonably  good  detection 
capability  out  to  80  km  in  range.  Stratospheric  returns  again  are  the 
highest  of  any  region;  however,  comparison  with  measured  stratospheric 
backscatter  data  indicates  that  the  dirty  atmosphere  model  is  not  very 
realistic . 

It  is  apparent  that  none  of  the  Ildars  are  capable  of  detecting 
wakes  in  the  clean  tropopause  and  free  tropospheric  regions.  This 
result  is  unfortunate  since  this  area  is  the  prime  location  for  cruising 
aircraft,  Also  evident  in  these  results  is  that  it  takes  a  major  volca¬ 
nic  eruption  to  provide  sufficient  signal  level  for  wake  detection  in 
the  stratosphere.  Obviously,  this  is  unacceptable.  One  interesting 
observation  is  that  while  there  is  n  good  correlation  between  backseat- 
ter  coefficient  and  extinction  coefficient  (they  both  increase 
together),  areas  with  higher  aerosol  backscatter  coefficient  showed  high 
signal  leve).<i  despite  the  Increase  in  extinction,  Thus,  backscatter 
coefficient  appears  to  be  the  veal  factor  for  determing  lidar  perform¬ 
ance  . 


4-37 


Surprisingly,  the  9.115  pm  lidar  performed  better  than  the  short 
wavelength  lidars  in  all  cases.  This  result  is  surprising  since  the 
aerosol  backscatter  coefficients  are  higher  at  the  shorter  wavelengths. 
There  may  be  several  explanations  for  this  observation.  First,  since 
shot  noise  power  is  inversely  proportional  to  wavelength,  shot  noise  is 
roughly  9  times  higher  at  1.069  Atm  than  at  9.115  pm.  Signal • to-noise 
ratio  is  thus  degraded  accordingly. 

Second,  loss  in  spatial  coherence  due  to  refractive  turbulence  at 
long  ranges  is  significant  at  1.069  pm  and  even  at  2.091  pm.  From 
Figure  2-9,  the  transverse  coherence  length  at  80  km  range  is  only  about 
0.13  m  for  the  1.069  pm  case.  Transverse  coherence  length  at  9.115  pm 
is  1.7  m  at  80  km  range,  still  much  larger  than  the  receiver  diameter. 
The  result  is  that  signal • to-noise  ratio  at  80  km  range  is  reduced  by  a 
factor  of  about  3  at  1.069  pm  and  by  a  factor  of  only  1.02  at  9.115  pm. 
These  calculations  assumed  a  constant  value  of  Cj  of  10'w  m'2'3,  which 
should  be  fairly  conservative. 

Third,  for  an  unfocused  system  the  wavefront  mismatch  at  short 
ranges  is  very  degrading  for  the  short,  wavelength  lidars  due  to  the  fact 
that  even  slight  curvature  of  the  received  wavefront  causes  destructive 
interference  with  the  plane  wave  local  oscillator  beam  at  the  edges  of 
the  detoctor.  The  term 

'nD*V 
.  4  \  A*  J 

in  the  signal  reduction  factor  <SRF)  describes  this  effect.  At  the 
relatively  short  range  of  20  km,  this  term  takes  on  a  value  of  11.0  at 
1.069  pm,  At  9.115  pm,  this  term  is  only  0.15.  The  degradation  in 
signal-to-noise  ratio  at  20  km  range  Is  thus  about  a  factor  12  at  1.069 
pm.  At  80  km  range,  this  degradation  is  down  to  only  1.7,  neglecting 
refractive  turbulence  losses. 
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Fourth,  the  total  attenuation  at  shorter  wavelengths  is  higher.  In 
addition  to  higher  aerosol  attenuation,  the  shorter  wavelengths  also 
receive  substantial  extinction  from  molecular  scattering.  At  80  km  hor¬ 
izontal  range  in  the  tropopause,  two  way  transmissions  at  1.064  pm  and 
9,115  pm  are  0.852  and  0.990,  respectively .  Fifth,  the  higher  molecular 
backscatter  at  shorter  wavelengths  also  degrades  signal-to-noise  ratio, 
although,  as  discussed  in  chapter  2,  this  effect  should  be  relatively 
small.  The  combination  of  these  five  explanations  appears  to  account 
for  the  relatively  poor  performance  of  the  short  wavelength  lidars  at 
all  ranges. 

It  seems  reasonable  to  conclude  that  high  altitude  wake  detection 
taking  into  account  only  naturally  occurring  aerosol  backscatter  is  not 
feasible  with  the  given  lidar  system  parameters;  however,  improvement  in 
signal-to-noise  ratio  may  be  made  by  using  higher  laser  pulse  energies, 
averaging  return  signals,  or  using  larger  receivers. 

Increasing  the  pulse  energy  of  the  transmitter  laser  will  also 
increase  signal-to-noise  ratio  proportionally.  For  this  analysis,  we 
assumed  a  pulse  energy  of  0.2  J.  By  increasing  pulse  energy  by  a  factor 
of  10  to  2  J,  signal -to-noisc  ratio  will  also  increase  by  a  factor  of 
10.  In  Figures  4-6  through  4-11,  the  resulting  minimum  signal  - to-noisc* 
ratio  would  be  -10  d&.  The  1.064  pm  lidar  still  would  not  be  capable  of 
any  appreciable  detection  range  for  all  but  the  dirtiest  of  atmospheres. 
The  performance  of  the  2.091  pm  lidar  would  be  much  improved,  but  long 
range  detection  in  the  tropopause  would  still  be  out  of  the  question. 

The  9.115  pm  lidar  would  perform  adequately  in  the  moderate  and  dirty 
atmospheres,  but  would  be  limited  in  the  clean  regions  of  the  atmo¬ 
sphere  . 

Receiver  area  is  also  directly  proportional  to  signal-to-noise 
ratio;  however,  refractive  turbulence  degradations  limit  the  maximum 
size  of  the  aperture.  Figure  2-9  shows  that  for  80  km  detection  range, 
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the  maximum  usable  receiver  diameters  are  0.13  m  at  1.064  pm,  0.3  m  at 
2.091  pm,  and  1.7  m  at  9.115  pm.  Obviously,  no  improvements  could  be 
made  in  the  1.064  and  2.091  pm  lidar  performances  at  ranges  beyond  80 
km.  Increasing  the  receiver  diameter  of  the  9.115  pm  lidar  to  1  m 
(probably  out  of  the  question  for  a  flyable  system!)  would  result  in  an 
increase  In  signal-to-noise  ratio  of  about  a  factor  of  11.  As  shown 
above,  this  would  still  leave  the  clean  arsaa  of  the  atmosphere  unde¬ 
tectable  . 

Probably  the  most  realistic  improvement  in  signal -to-noise  ratio 
could  come  from  averaging  the  return  signals  from  several  pulses  in  the 
same  scattering  volume.  The  ideal  improvement  in  signal- to-noise  ratio 
is  N1/a  for  N  independent  shots  averaged.  In  reality,  however,  the 
improvement  is  somewhat  less  due  to  "a  small  but  long-term  temporal 
autocorrelation  caused  by  time-varying  atmospheric  effects  (Menyuk  et 
al,,  1983:186)."  Thus,  averaging  10  shots  would  result  in  a  signal -to- 
noise  ratio  increase  slightly  lees  than  5  dB.  The  main  problem  with 
shot  averaging  is  the  scanning  limitation  that  it  imposes.  Since  the 
lidar  probe  beam  in  so  narrow,  in  order  to  scan  large  solid  angles,  a 
large  number  of  sample  pulses  would  be  required.  By  averaging  shots 
from  each  beam  angle,  the  required  pulse  rate  to  scan  a  reasonab1  solid 
angle  in  a  reasonable  amount  of  time  would  be,  well,  unreasonable. 

Even  with  the  signal-to-noise  ratio  improvements  discussed  above, 
aircraft  wake  detection  at  high  altitude  for  backscatter  from  atmo¬ 
spheric  aerosols  alone  does  not  appear  practical. 
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5.  Exhaust  Trails 

Particulates  in  the  exhaust  trail  of  an  aircraft  may  provide  an 
enhancement  in  backacatter  from  the  wake.  In  this  chapter,  a  first* look 
analysis  will  be  made  to  try  and  determine  the  magnitude  of  this 
enhancement.  Two  types  of  exhaust  clouds  will  be  studied:  1)  soot 
particles  from  engine  combustion,  and  2)  ice  particles  from  condensed 
water  vapor  in  the  exhaust.  Since  soot  is  a  by-product  of  combustion, 
it  should  always  be  present  at  aome  level  in  the  WAke.  During  certain 
atmospheric  conditions,  water  vapor  from  engine  exhaust  condenses  and 
forms  highly  visible  contrails.  It  la  expected  that  backacatter  from 
contrail  particles  is  very  high  and  will  easily  facilitate  lidnr  woke 
detection.  Backscatter  from  soot,  however,  may  or  may  not  be  substan¬ 
tial. 

5  +  L.  Sjoot-  Emissions 

Soot  particles  are  the  primary  particulate  in  the  hot,  uncondensed 
exhaust  from  Jet  engines.  Since  soot  is  always  present  in  the  exhaust 
trail,  it  may  provide  enhancement  in  the  backscatter  coefficient  within 
the  aircraft  wake.  In  the  following  analysis,  a  brief  description  of 
soot  formation  and  particle  characteristics  will  be  given.  Following 
that,  Mie  scattering  calculations  will  be  made  on  several  particle  size 
distributions  from  published  Jet  engine  exhaust  measurements,  From 
these  scattering  calculations,  backscatter  levels  in  the  equilibrium 
wakes  of  four  aircraft  will  be  estimated  to  determine  the  magnitude  of 
enhancement  in  lidar  wake  detection  performance, 

5.1.1. — laimtiaa 

Soot  is  the  by-product  of  imperfect,  or  incomplete,  combustion. 

Some  of  the  carbon  atoms  that  are  liberated  from  fuel  during  combustion 
do  riot  combine  with  oxygen  atoms  to  form  exhaunt  gases  such  as  CO  and 
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C02.  These  free  carbon  atoms  coagulate  together  to  form  soot  particles. 
Immediately  after  combustion,  most  of  the  soot  particles  are  composed  of 
clumps  of  around  60  carbon  atoms.  These  Initial  soot  particles  have  a 
particle  size  of  about  0.02  (Wander,  1990).  As  the  exhaust  flows  out 
of  the  engine  and  into  the  atmosphere,  soot  particles  collide  and  stick 
together  forming  larger  particles.  This  process  of  agglomeration  con* 
tlnues  as  the  exhaust  cloud  diffuses  into  the  atmosphere.  Unfortunately 
for  us,  these  particles  are  not  spherical  but  tend  toward  long  chains  of 
connected  smaller  particles  (Iskander  et  al.,  1989;  and  Wander,  1990). 

In  steady  state,  the  distribution  of  soot  particles  sizes  appears 
to  approach  a  lognormal  distribution  (Rosen  and  Greegor,  1974:244). 

Three  published  soot  size  distributions  are  listed  in  Table  5*1.  The 
AFGL  model  is  for  atmospheric  soot  from  urban  sources,  and  the  NRC  model 
for  is  soot  following  nuclear  detonation  in  the  atmosphere.  The  model 
by  Rosen  and  Greegor  is  the  model  fit  to  measured  particle  size  distrib¬ 
utions  in  the  wake  of  an  F-104  aircraft  at  30,000  feet  altitude.  These 
models  will  be  used  for  comparison  of  backscatter  coefficients 
calculated  from  measured  Jet  engine  exhaust  soot  distributions. 

Table  5-1.  Three  Published  Soot  Particle  Size  Distribution  Models 


Model 

Distribution 

Mean 

Radius 

0‘m) 

Standard 

Deviation 

Relative 

Number 

Density 

(cm'3) 

AFGL'' 

Sum  of  2  Log¬ 
normal* 

0.03 

0.5 

■ram 

■ram 

0.999875 

0.000125 

NRC"’ 

0.1 

0,30 

1 

Rosen  &  Greeg¬ 
or*** 

Sum  of  2  Log¬ 
normal  s 

0.08 

0,3 

■nra 

■oora 

— 

* ( Fenn  et  al. ,  1985:18-  ) 

**(NRC,  1985) 

•‘‘Adapted  from  (Rosen  and  Greegor,  1974:244) 
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Several  references  ware  located  that  contained  measured  particle 
size  distribution!  from  Jet  engine  exhaust  (Mathews  et  al.,  1984;  John¬ 
son  ee  al.,  1985;  and  Spicer  et  al.,  1987).  These  size  distributions 
were  input  to  computer  and  a  Mie  scattering  program  was  used  to 
calculate  the  total  scattering  coefficient  and  phase  functions  for  each 
distribution.  Using  the  equations  in  chapter  2,  the  volume  backscatter 
coefficient  was  then  computed.  Even  though  soot  particles  are  not 
spherical,  Mie  scattering  calculations  were  used  to  give  a  rough  esti¬ 
mate  of  the  magnitude  and  range  of  backscatter  coefficients  in  Jet 
engine  exhaust.  These  results  should  be  accurate  to  within  at  least  an 
order  of  magnitude,  if  not  better.  For  comparison,  backscatter  coeffi¬ 
cients  for  the  theoretical  soot  size  distributions  in  Table  5-1  were 
also  computed  using  the  Mie  code.  The  comparison  of  results  will  be 
discussed  later  in  this  section. 

The  measured  particle  size  distributions  are  listed  in  the  Appen¬ 
dix.  Table  5-2  lists  the  engine  type,  thrust  setting,  approximate  plume 
diameter  at  measurement  point,  and  total  particle  concentration  for  each 
data  set.  Some  extrapolations  of  the  size  distributions  were  necessary 
since  most  of  the  measurements  did  not  cover  a  wide  enough  range  of 
particle  sizes.  The  refractive  indices  of  soot  at  the  three  wavelengths 
of  interest  were  taken  from  Fenn  et  al.  (1985:18-17)  and  are  listed  in 
Table  5-3.  For  the  TF30-P1  and  J52-P3  data  sets,  the  particle  size 
distribution  data  was  given  in  mass  density  per  unit  volume.  For  the 
Mie  scattering  calculations,  these  data  were  converted  to  number  density 
using  the  equation 


n(r) 


3p  c 

Bn  Ar’p, 


(5-1) 
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where 


n(r)  -  particle  number  density  (cm'Vm'1) 

r  -  particle  radius  (cm) 
p,  -  particle  cloud  density  (g/cm3) 
p,-  soot  density  -  1,5  g/cm3 
A r  -  ra-r,  (cm) 

Table  5-2.  Jet  Engine  Exhaust  Data  Sets  Used  in  Mie  Scattering  Calcula¬ 
tions 


Engine  TypetMf) 

Thrust 

(lbs) 

Plume 

Diameter 

(m) 

Total  Particle 
Concentration 
(cm‘3) 

TF30-P11 

9,500" 

6 

9.78x108 

TF30-P11 

9,500" 

8" 

4.74x10* 

J52-P33 

8,500" 

6* 

2.35x10* 

J52  -P32 

OO 

O 

O 

» 

8* 

3.76x10* 

TF30-P1033 

518 

1.2* 

4 . 92x10* 

TF30-P1033 

2,934 

1,2* 

5.96x10* 

TF30-P1033 

7,317 

1,2* 

5.18x10* 

TF30-P1033 

9,785 

1,2" 

3.98x10* 

TF30-P1093 

563 

1.2* 

2.17x10* 

TF30-P1093 

3,200 

1,2" 

3.68x10* 

TF30-P1093 

7,998 

1,2" 

3.44x10* 

TF30-P1093 

10,633 

1,2* 

2.82x10* 

TF41-A23 

460 

1* 

3.31x10* 

TF41-A23 

3,818 

1* 

3.76x10* 

TF41-A23 

9,554 

1" 

3.32x10* 

TF41-A23 

12,937 

1" 

4.26x10* 

"Estimate 

^Mathews  et  al  .  ,  1984) 

J(Johnaon  et  al .  ,  1985) 

-  turbojet 


For  each  exhaust  data  set  above,  the  backscatter  coefficient  at 
each  of  the  three  infrared  wavelengths  was  computed  and  plotted  versus 
engine  thrust  divided  by  plume  area,  This  normalization  of  engine 
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Table  5-3.  Refractive  Index  of  Soot  at  Three  Infrared  Wavelengths 


Wavelength 

(^m) 

Refractive  Index 

1.064 

1 . 750-JO .440 

2.091 

1.805-jO. 495 

9.115 

2. 176-jO. 700 

thrust  allows  comparison  of  backscattar  coefficients  for  all  engines 
since  it  is  assumed  that  the  mass  of  soot  in  the  exhaust  plume  is  pro¬ 
portional  to  engine  thrust  and  is  inversely  proportional  to  plume  area. 
The  actual  relationship  between  engine  power  setting  and  backscatter 
coefficient  is  probably  complex  since  backscatter  is  very  dependent  upon 
the  particle  size  distribution.  Measured  data  seems  to  show  that  mean 
particle  size  increases  with  increasing  thrust,  while  the  total  number 
of  particles  does  not  change  appreciably  (Spicer  et  al.,  1987:60).  The 
data  is  shown  in  Figures  5-l.a-c.  These  graphs  show  that  backscatter 
coefficient  can  easily  span  two  orders  of  magnitude  between  clean  and 
dirty  engines.  Also  shown  on  the  graphs  are  power  law  best  line  fits  to 
data.  The  power  law  fit  has  no  known  theoretical  basis;  it  was  simply 
chosen  for  a  straight  line  fit  on  log-log  axes.  The  best  line  fits 
appear  Co  be  representative  of  moderately  clean  engines.  Newer,  clean 
burning  turbofan  engines  will  likely  fall  below  these  lines. 

From  chapter  2,  we  know  that  the  exhaust  products  becomo  wrapped  up 
in  the  recirculation  cell  surrounding  the  wake  vortices.  Since  no  mass 
is  transferred  across  the  recirculation  cell  boundary,  we  may  assume 
that  all  the  exhaust  particles  are  entrained  within  it.  In  addition,  we 
will  assume  that  these  particles  are  uniformly  distributed  within  the 
recirculation  cell.  This  assumption  seems  good  due  to  the  highly  turbu¬ 
lent  air  within  the  recirculation  cell  itself.  The  backscatter  coeffi¬ 
cient  in  the  wake  may  be  estimated  from  the  graphs  above  by  finding  the 
total  thrust  per  unit  area  in  r.he  recirculation  cell.  The  ratios  of 


5-5 


Backscatter  Coefficient 


2.091  /im 


Backscatter  Coefficient 


Normalized  Thrust  (lbs/m2) 


Figure  5*1. c.  Exhaust  Backscatter  Coef ficientes  at  9,115  versus  Nor 

malized  Thrust 


total  maximum  thrust  to  recirculation  cell  area  for  the  four  test  air¬ 
craft  are  listed  in  Table  5-4.  Table  5-5  lists  the  estimated  exhaust, 
backscatter  coeffients  in  the  wakes  of  these  aircraft  taken  from  the 
best  line  fits  to  the  data  in  Figures  5-1. a-c. 

Table  5-4.  Ratio  of  Maximum  Thrust  to  Recirculation  Cell  Area  for  Four 

Aircraft 


Aircraft 

Thrust/Area 

(lbs/m*) 

Large  Transport 

20 

Medium  Transport 

17 

Large  Fighter 

177 

Small  Fighter 

89 

Table  5-5.  Mean  Exhaust  Backscatter  Coefficients  in  the  Wakes  of  Four 

Aircraft 


Aircraft 

£.(*) . 

1. .  064  pm 
(m-ier*1) 

0.(ff)  . 

2,091  pm 
(m"lsr"l> 

9.115  pm 
(m^sr-1) 

Large 

Transport 

5xl0*7 

9x10*0 

2.2xlO-» 

Medium 

Transport 

4xl0"7 

7xl0-» 

1.8xl0-» 

Large 

Fighter 

1.7xl0-B 

1 .  Oxl'J-o 

2.0x10-' 

Small 

Fighter 

1.2xl0-« 

5xlO"7 

1 . 0xl0'e 

For  comparison,  the  three  model  size  distributions  in  Table  5-1 
were  also  run  through  the  Mie  scattering  program.  Since  these  models 

1 

represent  relatively  steady  state  soot  size  distributions,  they  should 

be  Indicative  of  the  soot  size  distributions  in  the  middle  and  late  wake  j 

stages.  The  total  number  densities  were  scaled  for  clean  and  dirty  : 

engines  based  on  total  soot,  mass  per  unit  volume,  which  am  readily  j 

available  parameters  of  jet  engine  soot  emissions.  Newer  high-bypass 

engines  have  soot  emissions  of  less  than  2  tng/in3  at  the  exit  nozzle,  ; 

while  some  older  engines  may  have  emissions  as  high  as  20-30  mg/m3 
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(Pitz,  1987:3).  Table  5-6  lists  the  calculated  backscatter  coefficients 
for  the  three  models  and  for  clean  (1  mg/m3  soot)  and  dirty  (30  mg/in* 
soot)  engines. 

Table  5-6.  Modeled  Soot  Backscatter  Coefficients  for  Clean  and  Dirty 

Jet  Engines 


Soot  Model 

Engine 

0.OO, 
1.064  pm 
(m^sr"1) 

0.00, 
2.091  pm 
(m-'sr"1) 

0.(ff)  , 
9.115  pm 
(m-’sr"1) 

AFGL 

Clean 

1.55x10'® 

9 . 17x10*’ 

2.76x10-’ 

II 

Dirty 

4.66x10-* 

2.75x10-* 

8.28x10-® 

NRC 

Clean 

1.52x10-® 

1.58x10-® 

2.65x10*’ 

II 

Dirty 

4.57x10'* 

4.73x10** 

7.94x10-® 

■mm 

■MB  WAJMi 

Clean 

1.38x10-* 

2.91x10-® 

3.34x10-’ 

li 

Dirty 

4.14x10-* 

8.74x10-* 

1.00x10-* 

The  data  above  shows  surprisingly  little  variance  from  model  to 
model.  This  result  lends  confidence  to  the  backscatter  calculations  for 
soot  since  backscatter  in  the  infrared  seems  to  be  only  mildly  sensitive 
to  changes  in  the  size  distributions.  This  insensitivity  is  probably 
due  to  the  small  particle  sizes  compared  to  the  infrared  wavelengths, 
used. 

In  order  to  compare  these  results  with  calculations  based  on  mea¬ 
sured  engine  particle  distributions,  some  assumptions  must  be  made  about 
engine  thrust  and  nozzle  diameter.  If  we  assume  engine  thrusts  from 
10,000  to  20,000  lbs  and  aim  nozzle  diameter,  then  the  normalized 
thrust  to  plume  areas  would  be  from  12,700  lbs/m2  to  25,400  lbs/m?. 

From  Figures  5-1. a-c,  the  backscatter  coefficients  for  modeled  data  are 
generally  lower  than  for  measured  data,  particularly  at  1.064  and  2,091 
pm.  If  the  normalized  thrust  for  the  modeled  size  distributions  was 
between  100-1000  lbs/m2,  the  data  would  correlate  well  with  measured 
data.  Only  measured  data  will  be  used  in  the  performance  analysis. 
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The  wake  backscatter  coefficients  in  Table  5*5  represent  the 
expected  levels  for  moderately  clean  engines.  These  values  may  be  com¬ 
pared  to  the  backscatter  coefficients  for  the  atmospheric  aerosol  models 
developed  in  chapter  4.  Except  for  boundary  layer  aerosols  and  extreme 
volcanic  aerosols,  soot  backscatter  in  the  four  aircraft  wakes  ero 
higher  than  the  ambient  aerosol  backscatter  levels. 

To  estimate  the  Improvement  in  Ildar  performance  due  to  the  pres¬ 
ence  of  exhaust  soot,  algnal-to-noise  ratio  was  computed  for  a  Ildar 
looking  horizontally  through  the  clean  atmosphere  tropopause  at  9  km 
altitude.  The  backscatter  coefficient  was  taken  to  be  the  sum  of  the 
ambient  aerosol  backscatter  coefficient  and  the  exhaust  backscatter 
coefficient.  Since  Ildar  performance  was  shown  to  be  worst  in  the  tro¬ 
popause  (chapter  4),  thia  eimulation  should  indicato  the  minimum  per¬ 
formance  assuming  enhancement  by  exhaust  soot.  Data  from  these 
simulations  are  shown  in  Figures  5*2. a-c.  These  simulations  do  not  take 
into  account  the  reduction  in  signal - to-noise  ratio  within  a  single 
range  gate  for  pulse  lengths  that  are  much  wider  than  the  wake. 

As  the  data  in  Figures  5-2. a-c  show,  the  enhancement  in  slgnnl-to- 
noise  ratio  due  to  exhaust  soot  is  very  significant  in  very  clean  atmo¬ 
sphere.  Interestingly,  the  fighter  aircraft  had  the  highest  enhancement 
due  to  the  much  higher  thrust/wake  area  ratios,  These  simulations  show 
detectability  of  fighter  aircraft  wakes  in  excess  of  80  km  In  the  tropo- 
pnuse  at  all  wavelengths.  Maximum  detection  range  for  the  transport 
aircraft  varies  from  50-75  km.  If  the  exhaust  backscatter  coefficients 
used  here  are  reduced  by  a  factor  of  10,  the  corresponding  signal -to- 
noise  ratio  will  also  be  reduced  by  a  factor  of  10.  This  reduction 
would  eliminate  any  wake  detection  by  the  1.064  pm  as  well  as  transport 
wakes  by  the  2,091  pm  lidar.  Maximum  detection  range  of  fighter  air- 
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Signal— to— Noise  Ratio 


Range  (km) 


Figure  5*2, a,  Signal*to*Noise  Ratio  in  the  Soot  Enhanced  Wakes  of  Four 
Aircraft  at  Ranges  from  0-80  km.  Wavelength  is  1,06^*  and  lidar  is 
at  9  km  altitude  looking  horizontally. 
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Signal— to- Noise  Ratio 


Figure  5 - 2  . b 
Aircraft  at 


100 


10d 


Laras  Transport 
Mtalum  Transport 
Largs  Flghtsr 
Small  Flghtsr 


0.1 


Range  (km) 


Signal-to-Noine  Ratio  in  the  Soot  Enhanced  Wakes  of  Four 
Ranges  from  0>80  km.  Wavelength  is  2.091  ^m,  and  lidar  is 
at  9  km  altitude  looking  horizontally. 
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Signal— to— Noise  Ratio 


Figure  5-2. c.  Signal- to-Noise  Ratio  in  the  Soot  Enhanced  Wakes  of  Four 
Aircraft  at  Ranges  from  0-80  km,  Wavelength  is  9,115  ,  and  lidar  is 

at  5  km  altitude  looking  horizontally. 


craft  wakes  would  fall  to  50-80  km  at  2,091  /uni  and  40-60  km  at  9,115  /um , 
This  is  still  an  Improvement  over  detection  performance  for  atmospheric 
aerosol  backscatter  only. 

An  additional  benefit  of  the  enhancement  by  exhaust  particles  is 
the  fact  that  velocity  width  of  the  return  signal  will  not  be  degraded 
as  much  by  longer  pulse  lengths.  Figures  3-5. a-d  showed  how  this  wake 
signature  was  reduced  at  longer  pulse  lengths;  however,  these  graphs 
were  produced  assuming  backscatter  was  homogeneous  inside  end  outside 
the  wake.  For  exhaust  enhanced  wake  signals,  the  highly  turbulent  air 
within  the  recirculation  cell  will  contribute  more  to  the  signal  than 
the  air  outside  this  region.  The  resulting  velocity  width  of  the  return 
signal  should  be  higher  and  the  wake  should  be  even  easier  to  detect 
than  for  the  ambient  aerosol  case. 

UL _ Contrails 

In  contrast  to  soot,  contrails  are  not  always  presont  in  aircraft 
wakes;  however,  due  to  their  high  viaiblliy,  it  is  assumed  that  the 
enhancement  in  backscatter  from  aircraft  contrails  will  be  very  signifi¬ 
cant.  In  this  subsection,  aircraft  contrail  formation  will  be  briefly 
described,  estimations  will  be  made  of  the  backscatter  levels  from 
contrails,  and  corresponding  detection  performance  will  be  annlyz.u. 

5  ,-2.1-, _ Esiinuillan 

A  technical  reference  put  out  by  the  Air  Force  Air  Weather  Service 
gives  a  good  overview  of  contrail  formation  (Air  Weather  Service, 

1981:1) : 

Condensation  trails  (contrails)  are  elongated  tubular- 
shaped  clouds  composed  of  water  droplets  or  ice  crystals  which 
form  behind  aircraft  when  the  wake  becomes  supersaturated  with 
respect  to  water. 
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The  most  important  type  of  condensation  trail  arises  when 
the  water  vapor  in  the  exhaust  gas  mixes  with  and  saturates  the 

air  in  the  wake  of  the  aircraft  (Appleman,  1953),  Combustion 

of  the  hydrocarbon  fuels  used  in  aircraft  •  both  propeller  and 
jet  •  injects  both  water  vapor  and  heat  into  the  wake.  The 
added  moisture  raises  the  relative  humidity  in  the  wake,  while 
the  added  heat  lowers  it.  Whether  or  not  the  wake  will  reach 
saturation  depends  on  the  ratio  of  water  vapor  to  heat  in  the 

exhaust  gas  and  on  the  initial  pressure,  temperature,  and  rela* 

tive  humidity  of  the  environment. 

Figure  5-3,  taken  from  the  same  source,  shows  a  graph  of  the  rela¬ 
tive  humidity,  pressure,  and  temperature  required  for  contrail  formation 
from  Jet  aircraft.  The  graph  is  for  theoretical  behavior;  actual 
contrail  formation  varies  from  this  predicted  behavior  slightly  (Air 
Weather  Service,  1981:10-11).  The  diagonal  lines  across  the  graph  rep¬ 
resent  the  minimum  relative  humidirv  -squired  for  contrail  formation. 
Regions  to  the  left  of  the  0%  relative  humidity  line  are  where  contrails 
should  always  form,  and  regions  to  the  right  of  the  100%  relative  humid¬ 
ity  line  are  where  contrails  should  never  form.  Note  that  at  0%  rela¬ 
tive  humidity,  contrails  may  still  form,  but  it  takes  a  lower 
temperature  than  for  100%  relative  humidity  (at  constant  pressure), 

Also  of  note  is  the  fact  that  engine  power  setting  is  nearly  independent 
of  contrail  formation.  Only  contrail  intensity  varies  with  engine  power 
(Air  Weather  Service.  1981:8). 

In  a  brief  literature  search  on  contrails,  no  particle  size  distri¬ 
bution  data  was  located.  Additionally,  no  sources  were  found  that  indi¬ 
cated  when  contails  contained  water  droplets  and  when  they  contained  ice 
crystals.  Contrails  are  likely  to  be  composed  of  mostly  water  droplets 
immediately  after  formation.  After  some  period  of  time,  these  water 
droplets  may  form  into  ice  clouds  similar  in  form  to  cirrus  clouds  if 
conditions  are  right.  A  more  in-depth  search  is  obviously  needed, 
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Figure  5'3,  Jet  Aircraft  Contrail 
Humidity,  Pressure,  and  Temperat 


Since  no  size  distribution  data  were  located,  some  assumptions  were 
made  on  size  distribution  and  number  density  in  order  no  calculate  back- 
scatter  coefficients.  It  was  assumed  that  contrails  would  be  primarily 
composed  of  ice  crystals.  This  assumption  should  not  have  a  great 
impact  on  final  results  since  the  backscatter  coefficients  from  water 
droplet  contrails  should  not  be  much  different  as  the  total  water  con¬ 
tent  would  be  the  same.  Thus,  cirrus  cloud  particle  size  distributions 
were  used.  The  size  distributions  for  four  types  of  cirrus  clouds  were 
used:  1)  cirrostratua ,  2)  cirrus  uncinus,  3)  warm  cirrus,  and  4)  cold 
cirrus  (Takano  and  Liou,  1989:6).  The  appendix  lists  the  size  distribu¬ 
tion  data  for  each  cirrus  model. 

Although  ice  particles  are  not  spherical,  these  size  distributions 
were  run  through  the  Mie  scattering  program  to  calculate  total  scatter¬ 
ing  coefficients  and  phase  functions.  The  complex  refractive  indices  of 
ice  at  the  three  wavelengths  of  interest  are  listed  in  Table  5-7 
(Warren,  1984).  For  comparison,  Che  scattering  coefficients  derived 
from  Mie  scattering  calculations  are  listed  in  Table  5-8  along  with  the 
scattering  coefficients  at  1  and  2  n m  derived  by  Takano  and  Liou 
(1989:9)  for  hexagonal  ice  crystals.  Unfortunately,  no  9  jjm  data  was 
located.  Finally,  the  backscatter  coefficients  from  Mie  scattering 
results  are  shown  in  Table  5-9, 

The  minimum  cutoff  wavelength  used  in  the  Mie  calculations  was  0.01 
*im;  however,  no  mention  was  made  in  the  source  article  about  minimum 
cirrus  particle  size.  To  test  the  sensitivity  of  the  results  to  larger 
minimum  particle  sizes,  the  Mia  code  was  run  again  on  the  cirrus  uncinus 
model  for  minimum  radii  of  0.01  0.1  jim,  1.0  jim,  and  10.0  jim.  At  a 

minimum  size  of  10.0  jim,  the  1.064  pm  wavelength  backscatter  coefficient 
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dropped  by  «  factor  of  1.7.  The  9.115  /iin  wavelength  was  totally  insen¬ 


sitive  to  these  changes.  The  overall  results,  therefore,  should  not  he 
affected  very  much  by  this  assumption  on  minimum  particle  sir.e. 

Table  5-7.  Refractive  Index  of  Ice  at  Three  Infrared  Wavelengths 


Wavelength 

(Mm) 

Refractive  Index 

1 .064 

1.301-J1. 96x10-9 

2.091 

1 . 270-J9 . 0x10*' 

9,115 

1.265-J4. 45x10*2 

Table  5-8.  Comparison  of  the  Ice  Crystal  Scattering  Coefficients 
derived  from  Mie  Scattering  and  by  Takano  and  Liou  for  Hexagonal 

crystals 


Cirrus  Cloud  Model 

Wavelength 

(Mm) 

0,. 

Mlo  Scattering 
(km*1) 

I'aknno  &  U  mi 
(km*1) 

Cirrostratus 

1 

0,3969 

0.3863 

II 

2 

0,3504 

0.3550 

Cirrus  Uncinus 

1 

5.792 

2.601 

II 

2 

3.646 

2,084 

Warm  Cirrus 

1 

0,8215 

0,6521 

II 

2 

0.6312 

0.5973 

Cold  Cirrus 

1 

0.0963 

0.1662 

li 

2 

0.0848 

0. 1584 

Table  5-9.  Backscntter  Coefficients  for  Four  Cirrus  Cloud  Models 

Derived  from  Mie  Scattering 


Cirrus  Model 

UMO. 

1 . 064  /ini 
(nr'sr-1) 

0.(0 

2.091  MU' 
(nr'ar*1 ) 

0.(0 

9.115  /an 
( in* 1  sr1 ) 

Cirrostratus 

l .  2x10*° 

4.0xl0-J 

1.5x10-8 

Cirrus  Uncinus 

2 . 9xl0's 

1 .  Ixl0-B 

2.5x10*’ 

Warm  Cirrus 

4 . 8x1 0  B 

3,6x10-’ 

5.0x10" 

Cold  Cirrus 

3 , 7x10"’ 

2.9xlO-« 

1.2x10" 

The  results  in  Table  5-9  carry  a  number  of  assumptions  with  them. 
They  assume  spherical  particles  that  are  randomly  oriented  and  spaced. 
In  addition,  they  assume  the  particles  are  isotropic  and  that  the  hack- 
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Table  5-10,  Measured  Cirrus  Cloud  Backscatter  Coefficients  at  10,6  nm 


Koenig  and  Alejandro,  1990 
2Hal 1  et  al. ,  1988 
3Gross  et  al, ,  1984 

scattered  light  is  of  the  same  polarization  as  the  incident  l<ght. 
Despite  all  of  these  assumptions,  the  results  above,  at  least  at  9.115 
/am,  fare  reasonably  veil  with  measured  data  at  10.6  /am.  Table  5-10 
shows  cirrus  cloud  backscatter  coefficients  measured  at  10.6  /um  versus 
cloud  altitude  for  a  number  of  different  cirrus  clouds,  The  data  is 
combined  from  several  sources,  though  much  of  it  was  taken  during  a 
Joint  U.S.  Air  Force/U.K.  Royal  Signals  and  Radar  Establishment  measure¬ 
ment  program  above  Ascension  Island  (Koenig  and  Alejandro,  1990;  Hall  et 
al.,  1988;  and  Gross  et  al.,  1984).  The  wide  range  of  cirrus 
backscatter  coefficients  Indicates  the  variability  in  cirrus  cloud 
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thickness  and  density.  Though  this  data  is  by  no  means  a  statistically 
complete  data  9et,  it  shows  the  magnitude  and  range  of  cirrus  backscat- 
ter  coefficients  at  10.6  pm.  The  geometric  mean  backscatter  coefficient 
from  this  data  is  6.9xl0‘s  m^sr*1,  and  the  geometric  standard  deviation 
is  12.6.  The  mean  value  from  this  data  la  about  a  factor  of  5  lower 
than  the  calculated  data,  wh<ch  has  a  geometric  mean  of  3.9x10  4  m^sr'1 
at  9.115  pm,  but  this  is  well  within  one  standard  deviation.  In  addi* 
tion,  based  on  the  trend  indicated  in  calculated  results,  cirrus  back¬ 
scatter  coefficient  decreases  with  increasing  wavelength  in  the  infrared 
region. 

The  backscatter  from  within  a  contrail  is  expected  to  be  highM* 
than  backscatter  from  naturally  occurring  cirrus  clouds.  The  basis  of 
this  assumption  is  the  dense  appearance  of  contrails  compared  to  nearby 
cirrus  clouds.  Some  Ildar  data  exists  to  verify  this  assumption.  A 
measurement  by  Koenig  using  a  0.532  pm  lidar  shows  a  strong  contrail 
return  at  about  10.8  km  altitude  above  a  thin  cirrus  layer  at  10.0  km, 
The  contrail  signal  was  roughly  an  order  of  magnitude  stronger  than  the 
cirrus  signal  and  was  probed  about  15  minutes  after  the  contrail  was 
formed  (Koenig,  1990).  The  contrail  thickness  was  less  than  150  m  ver¬ 
tically,  This  data  was  the  only  lidar  contrail  measurement  that  could 
be  found.  Based  on  this  evidence,  a  conservative  estimate  of  contrail 
backscatter  coefficient  would  be  10  times  higher  than  the  geometric  mean 
cirrus  backscatter  coefficient  from  the  four  cirrus  models  above.  Based 
on  this  assumption,  Table  5-11  gives  the  estimated  contrail  backscatter 
coefficients  at  the  three  infrared  wavelengths  of  interest. 

Table  5-11.  Estimated  Contrail  Backscatter  Coefficients 


Wavelength 

(nr’sr*1) 

1,064  pm 

2,8xlO-5 

2.091  pm 

2.6x10-® 

9.115  pm 

3 . 9xl0"7 
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As  for  the  soot  analysis,  signal-to-noise  ratio  for  each  lidar  was 
computed  for  a  lidar  at  9  km  altitude  looking  horizontally.  The  back- 
scatter  coefficient  at  this  altitude  was  taken  as  the  sum  of  the  atmo¬ 
spheric  aerosol  backscatter  coefficient  and  the  estimated  contrail 
backscatter  coefficient  from  Table  5-9,  The  results  are  plotted  in 
Figure  5*4  for  contrails  from  0-200  km  in  range. 

As  expected,  signal-to-noise  ratio  from  contrails  apt  ars  to  be 
extremely  high.  Figure  5-4  implies  contrail  detectability  well  beyond 
the  200  km  limit  of  the  graph.  Again,  results  indicate  that  the  9.115 
im  lidar  performed  best,  probably  due  to  the  degrading  refractive  turbu¬ 
lence  effects  at  the  shorter  wavelengths.  The  1.064  pm  and  2.091  n 
lidars  had  nearly  equal  performance  in  this  simulation.  Even  if  the 
contrail  backscatter  values  are  a  factor  of  10  too  high,  all  of  the 
Ildars  would  still  be  capable  of  detecting  the  contrail  enhanced  wakes 
at  horizontal  ranges  beyond  100  km. 


Range  (km) 


Figure  5-4,  Signal-to-Noise  Ratio  for  Contrail  Enhanced  Wakes  at  Range 
from  0-200  Tcni,  Lidar  is  at  9  km  altitude  looking  horizontally. 


6.  Conclusions 


This  thesis  has  analyzed  the  performance  of  three  infrared  Doppler 
lidar  systems  for  detecting  aircraft  wake  vortices  at  high  altitudes. 
Detection  performance  assuming  backscatter  from  atmospheric  aerosols 
alone  was  analyzed.  The  enhancements  in  backscatter  and  resulting  wake 
detectability  resulting  from  exhaust  particulates  were  then  analyzed.  A 
short  summary  of  these  results  and  recommendations  for  further  research 
are  outlined  below. 

6.1.  Summary  of  Results 

Probably  the  easiest  wake  signature  to  detect  is  the  increase  in 
spectral  width  of  the  return  signal  from  within  turbulent  wake  regions, 
An  analysis  of  the  spread  in  spectrum  width  for  four  aircraft  wakes 
revealed  that  a  minimum  signal* to-noite  ratio  of  about  0  dB  was  neces¬ 
sary  to  the  detect  the  wakes  in  all  cases.  This  minimum  value  was  used 
in  the  subsequent  analyses. 

Atmospheric  aerosols  alone  are  not  sufficient  for  detecting  air¬ 
craft  wakes  at  ranges  up  to  80  km  in  clean  atmospheres.  The  improve¬ 
ments  in  signal- to-noisa  ratio  discussed  at  the  end  of  chapter  4  are 
either  impractical  or  Insufficient  to  bolster  performance  in  very  clean 
regions  of  the  atmosphere.  Interestingly,  the  9.115  /jit.  lidar  performed 
best  in  simulations.  For  the  shorter  wavelengths,  the  relative  signal- 
to-noise  ratio  degradations  are  due  primarily  to  increased  shot  noise, 
wavefront  mismatch  at  short  ranges,  refractive  turbulence  at  long 
ranges,  and  higher  overall  extinction  from  both  aerosols  and  molecules. 

Enhancement  in  backscatter  due  to  engine  exhaust  trails  does  appear 
to  be  sufficient  for  wake  detection  beyond  80  km,  even  in  clean  atmo¬ 
sphere.  The  fighter  aircraft  ahowod  the  greatest  enhancement  in  wake 
backscatter  due  to  the  higher  thrust  to  wake  area  ratio  of  these 
aircraft.  For  the  soot  backscatter  coefficients  calculated  in  chapter 


6-1 


S ,  detection  of  fighter  aircraft  wakes  well  beyond  80  km  may  be  possible 
at  all  three  wavelengths.  There  la,  of  course,  some  uncertainty  in  the 
backscatter  estimates;  however,  even  for  backscatter  coefficients  that 
are  a  factor  of  ten  less  than  these  estimates,  detection  of  fighter 
aircraft  wakes  appear  possible  at  ranges  from  40-80  km  by  both  the  2.091 
jim  and  9.115  pm  Ildars. 

As  expected  the  estimated  Improvement  in  performance  due  to  the 
presence  of  contall  particles  was  very  significant.  All  Ildars  showed 
wake  detectability  for  all  four  aircraft  analyzed  at  ranges  beyond  200 
km  when  contrails  were  present.  Even  for  contrail  backscatter  values  a 
factor  of  10  lower  than  estimated,  wake  detection  was  still  possible 
beyond  100  km  in  range  for  all  cases. 

Currently,  there  is  a  great  deal  of  excitement  about  the  prospect 
of  coherent  lidar  using  solid  state  lasers  at  1-2  jim.  The  poor  signal- 
to-noiae  ratio  revealed  in  this  study  for  the  1.064  pm  and  2.091  pm 
Ildars,  however,  indicates  that  longer  wavelength  lidars,  such  as  C02, 
may  actually  perform  better  despite  the  higher  backscatter  coefficients 
at  1-2  pm. 

fi-J L RafiamiMndaLiana 

The  results  for  the  atmospheric  aerosol  and  contrail  enhanced  wakes 
seem  to  be  fairly  conclusive;  however,  the  results  for  the  soot  enhanced 
wakes,  while  appearing  positive,  are  not  completely  conclusive.  It 
appears  that  soot  emissions  will  provide  the  necessary  enhancement  in 
backscatter  for  long  range  detection,  but  this  analysis  made  a  lot  of 
assumptions.  Further  research  Into  the  scattering  properties  of  Jet 
engino  soot  are  needed  and  should  ideally  be  supported  by  lidar  measure¬ 
ments  . 
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Further  study  of  detection  methods  for  wake  return  signals  is  prob¬ 
ably  also  warranted.  Wake  signatures  may  exist  that  could  be  identified 
from  signals  below  0  dB,  Such  a  finding  would  do  much  to  improve 
detection  range. 
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Table  A-l.  Soot  Site  Dlatributlona  for  TP30-P1  Engine  Data 


Particle 

Teat  #1 

Teat  #3 

Radiua 

Number 

Number 

(pm) 

Denalty 

Denalty 

(cm^pm**) 

(cnf*pnf») 

0.014 


0.047 


0.080 


0.173 


0.342 


0.531 


0.837 


1.283 


2.333 


4.163 


1 . 3E+10 


1.58E+9 


1.58248 


4.10E+6 


1.65E45 


9.23E+3 


2.30E+3 


.0 


0.0 


6.26E49 


7.50E48 


9.30E47 


4.45E46 


1.57E45 


9.02E43 


2.26E43 


6.420 

0.592 

9.230 

0.352 

13.87 

0.0 

0.576 


0.160 


Table  A-2.  Soot  Site  Dietributione  for  J52-P3  Engine  Data 


Particla 

Teat  #7 

Teat  #6 

Radiua 

Number 

Number 

(pm) 

Density 

Denaity 

(cm‘5pm"i) 

(cm*3pnr*) 

0.007 

1.B6E48 

0.012 

1.63E47 

0,021 

2.34E47 

0,038 

0.0 

0.067 

6.96E45 

0.119 

1.B7E45 

0.211 

5.25E44 

0.375 

1.56E+4 

1.08E+8 


1.78E48 


6.11E+6 


1.63E46 


3.54E45 


1.63E44 


1.98E+4 


I 

I 


Table  A- 3 .  Soot  Size  Distributions  for  TF30-P103  Engine  Data 


Table  A-4.  Soot  Size  Distributions  for  TF30-P109  Engine  Data 


Tablo  A- 5.  Soot  Size  Distributions  for  TF41-A2  Engine  Data 


Particle 

Radius 

<^m) 

Idle 

Number 

Density 

(crn’tym'*) 

301  Power 
Number 
Density 
(crn'tym'*) 

75%  Power 
Number 
Density 
(cnftym**) 

lOOt  Power 
Number 
Density 
(cm‘tyra‘‘i) 

0,002 

3 . 1E+6 

0,0 

0.0 

0.0 

0.004 

4.5E+8 

0.0 

0.0 

3.3E+3 

0.007 

1.0E47 

0.0 

— 3E— 

7 . 8E+3 

0.012 

4.6E+6 

7.4E+7 

HEM 

0.0 

0.021 

9.3E+7 

1.0E+8 

0.0 

0.0 

0.038 

1.2E+7 

2.0E+7 

2 . 6E+6 

0.0 

0.067 

6.7E+5 

1.0E+7 

6.3E+7 

2.3E+7 

0.119 

8.8E+4 

9.2E+5 

2.7E+6 

1.9E+7 

0.3 

1.7E+5 

7.2E+5 

2.8E+6 

4.9E+6 

Table  A*6.  Ice  Crystal  Size  Distributions  for  4  Cirrus  Cloud  Models 


Particle 

Radius 

(/um) 

Cs* 

Number 

Density 

(cm'tym-1) 

Ci  Uneinus* 
Number 
Density 
(cm  **pnr') 

Warm  Ci* 
Number 
Density 
(cnrtynf*) 

Cold  Ci* 
Number 
Density 
(cm'Vnf  *) 

9 

1.7E-3 

1 .4E-3 

5.3E-3 

2.2E-3 

25 

2.5E-4 

3.1E-4 

1.8E-4 

1.9E-3 

65 

8.0E-6 

3.0E-5 

8.8E-6 

2.1E-7 

158 

3.9E-8 

') ,  8F.-6 

1.2E-6 

1.4E-8 

475 

0.0 

1.4E-6 

3.5E-8 

0.0 

"Cs  -  cirrostratus ,  Ci  •  cirrus 
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